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ABSTRACT 


A  wireless  sensor  network  (WSN)  consists  of  a  large  number  of  small  sensor 
nodes  that  are  densely  deployed  over  an  area  to  acquire  information  about  targets  of  in¬ 
terest.  These  sensor  nodes  collaborate  among  themselves  to  form  an  ad-hoc  network  and 
disseminate  the  collected  target  information  to  an  unmanned  aerial  vehicle  (UAV).  The 
objective  is  then  to  increase  the  data  rate  and  transmission  range  between  the  sensor 
nodes  and  the  UAV.  A  distributed  beamforming  approach  was  proposed  whereby  the 
sensor  nodes  are  grouped  into  clusters  and  their  transmission  are  coordinated  in  order  to 
form  a  distributed  antenna  array  that  directs  a  beam  towards  the  UAV. 

A  simulation  model  was  developed  and  implemented  in  MATLAB  programming 
language  to  study  the  effectiveness  of  beamforming  using  sensor  clusters  for  establishing 
a  communication  link  to  the  UAV.  Results  showed  that  the  antenna  main  lobe  remained 
relatively  unchanged  in  the  presence  of  position  errors  and  sensor  node  failures  or  when 
the  density  of  the  sensor  nodes  changed.  Additionally,  the  maximum  average  power  gain 
of  the  main  lobe  can  be  increased  by  increasing  the  density  of  the  sensor  cluster,  thereby 
increasing  the  transmission  range  between  the  sensor  clusters  and  the  UAV. 
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EXECUTIVE  SUMMARY 


A  wireless  sensor  network  (WSN)  consists  of  a  large  number  of  small  sensor 
nodes  that  are  densely  deployed  over  an  area  to  acquire  information  about  targets  of  in¬ 
terest.  They  can  be  deployed  for  remote  sensing  of  nuclear,  biological,  and  chemical 
weapons,  potential  terrorist  attack  detection  and  reconnaissance.  These  sensor  nodes  col¬ 
laborate  among  themselves  to  form  an  ad-hoc  network  and  disseminate  the  collected  tar¬ 
get  information  to  an  unmanned  aerial  vehicle  (UAV).  Advances  in  sensor  and  UAV 
technologies  result  in  higher  data  volume  transfer  and  longer  transmission  range  between 
the  sensor  nodes  and  the  UAV.  The  objective  is  then  to  increase  the  data  rate  and  trans¬ 
mission  range  between  the  sensor  nodes  and  the  UAV.  A  distributed  beamforming  ap¬ 
proach  was  proposed  whereby  the  sensor  nodes  are  grouped  into  clusters  and  their  trans¬ 
mission  are  coordinated  in  order  to  form  a  distributed  antenna  array  that  directs  a  beam 
towards  the  UAV. 

A  simulation  model  was  developed  and  implemented  in  the  MATLAB  program¬ 
ming  language  to  study  the  application  of  beamforming  using  sensor  clusters  for  estab¬ 
lishing  a  communication  link  to  the  UAV.  Different  sensor  node  densities  (i.e.,  number 
of  sensor  nodes  per  unit  area),  position  errors  and  sensor  node  failures  were  simulated  to 
investigate  their  effects  on  the  antenna  beam  generated  by  the  sensor  cluster.  Results 
showed  that  the  average  sidelobe  levels  increased  while  the  maximum  average  power 
gain  decreased  when  the  density  of  the  sensor  nodes  decreased.  However,  the  shape  of 
the  main  lobe  remained  relatively  unchanged.  Results  also  showed  that  the  presence  of 
position  errors  and  sensor  node  failures  reduced  the  maximum  average  power  gain  of  the 
antenna  beam  and  increased  its  sidelobe  levels,  but  the  shape  of  the  main  lobe  remained 
relatively  unchanged.  These  simulation  results  were  found  to  be  in  good  agreement  with 
the  theoretical  results  for  random  arrays. 

Since  the  shape  of  the  antenna  main  lobe  remained  relatively  unchanged  in  the 
presence  of  position  errors  and  sensor  node  failures  or  when  the  density  of  the  sensor 
nodes  changed,  it  can  be  steered  towards  the  UAV  in  an  adaptive  manner  as  the  UAV 
flew  over  the  sensor  field  independent  on  the  sensor  node  density,  position  errors  and 

xvii 


sensor  node  failures.  This  resulted  in  a  higher  signal-to-noise  ratio,  hence  a  higher  chan¬ 
nel  capacity.  Additionally,  the  maximum  average  power  gain  of  the  main  lobe  can  be  in¬ 
creased  by  increasing  the  density  of  the  sensor  cluster,  thereby  increasing  the  transmis¬ 
sion  range  between  the  sensor  clusters  and  the  UAV. 
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I.  INTRODUCTION 


Advances  in  wireless  communications  and  electronics  have  enabled  the  develop¬ 
ment  of  low-cost,  low-power  multifunctional  sensor  nodes  that  are  small  in  size  [1]. 
These  smart  disposable  microsensors  can  be  deployed  in  the  battlespace,  on  the  ground, 
in  the  air,  underwater,  in  vehicles,  inside  buildings  and  on  human  bodies  [2].  A  wireless 
sensor  network  (WSN)  consists  of  a  large  number  of  small  sensor  nodes  that  are  densely 
deployed  within  an  area  of  interest.  The  initial  exact  position  of  the  sensor  nodes  need 
not  be  predetermined  and  they  can  for  example  be  deployed  randomly  by  an  unmanned 
aerial  vehicle  (UAV)  as  shown  in  Fig.  1. 

z-axis 


Figure  1 .  Random  deployment  of  sensor  nodes  by  an  UAV. 


These  sensor  nodes  can  be  deployed  quickly  and  unattended  by  humans,  which 
make  them  very  attractive  for  military  applications.  They  can  be  deployed  for  remote 
sensing  of  nuclear,  biological,  and  chemical  weapons,  potential  terrorist  attack  detection 
and  reconnaissance  [1].  The  WSNs  provide  the  option  of  placing  sensor  nodes  very  close 
to  the  intended  targets  of  interest  without  putting  human  lives  at  risk. 
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Once  deployed,  these  multifunctional  sensor  nodes  have  the  ability  to  acquire  in¬ 
formation  about  targets  of  interest  and  collaborate  among  themselves  to  form  an  ad-hoc 
network  and  disseminate  the  collected  target  information  to  a  relay  node.  Figure  2  shows 
how  an  UAV  can  receive  the  data  from  these  sensor  nodes  in  a  wireless  manner  as  the 
UAV  flies  across  the  sensor  field.  The  sensor  nodes  first  detect  the  presence  of  the  UAV, 
which  approaches  from  an  elevation  angle  of  60  and  an  azimuth  angle  of  j>o  before  trans¬ 
mitting  the  data  to  the  UAV. 


Figure  2.  Wireless  data  transfer  from  the  sensor  nodes  to  the  UAV. 

A.  THESIS  OBJECTIVE 

The  operation  scenario  shown  in  Figures  1  and  2  is  adopted  in  this  thesis.  As  the 
sensor  technologies  improve,  more  sensor  components  can  be  incorporated  into  the  sen¬ 
sor  node,  which  results  in  more  sensing  data  being  collected  and  higher  amounts  of  data 
being  transmitted  to  the  UAV.  Additionally,  as  the  technology  improves,  the  UAVs  can 
fly  at  higher  altitudes  and  maneuver  more  efficiently.  The  objective  of  this  work  was  to 
increase  the  data  rate  and  transmission  range  between  the  sensor  nodes  and  the  UAV  in 
order  to  handle  the  higher  data  rates  and  longer  transmission  ranges  required  for  future 
military  applications. 
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The  ability  to  handle  the  increase  in  data  rate  and  transmission  range  between 
each  sensor  node  and  the  UAV  cannot  be  achieved  easily  because  these  sensor  nodes  are 
battery-operated  and  have  limited  transmission  power.  In  this  thesis,  a  distributed  beam¬ 
forming  approach  is  proposed  whereby  the  sensor  nodes  are  grouped  into  clusters  and 
their  transmissions  are  coordinated  in  order  to  form  a  distributed  antenna  array  that  di¬ 
rects  a  beam  towards  the  UAV  [6]. 

The  main  advantages  of  using  a  sensor  cluster  to  form  an  array-like  transmission 
to  the  UAV  are  increased  data  rate  and  transmission  range.  Since  the  arrays  can  form  a 
directed  beam  towards  the  UAV,  the  radiated  energy  is  aimed  in  the  direction  of  the 
UAV.  This  increases  the  signal-to-noise  ratio  at  the  UAV  and  results  in  increased  chan¬ 
nel  data  rate  as  given  by  Shannon’s  formula  [  3].  Additionally,  the  fact  that  a  sensor  clus¬ 
ter  is  a  type  of  multiple-input-multiple-output  (MIMO)  system  can  be  exploited  to 
achieve  further  data  rate  gains  [4].  The  antenna  gain  obtained  from  an  array  of  elements 
can  be  of  several  orders  higher  than  that  of  a  single  element.  This  additional  antenna  gain 
can  extend  the  transmission  range  between  the  sensor  nodes  and  the  UAV  for  a  fixed  sig- 
nal-to-noise  ratio  at  the  UAV.  Other  advantages  of  using  arrays  include  the  capability  to 
perform  spatial  filtering  by  steering  the  antenna  beam  towards  the  UAV,  thereby  sup¬ 
pressing  interference  signals  from  other  directions.  Additionally,  arrays  provide  the 
flexibility  to  control  the  antenna  gain  and  beamwidth  by  changing  the  number  of  array 
elements.  These  potential  benefits  make  the  sensor  cluster  a  suitable  candidate  to  meet 
the  objective  stated  earlier. 


B.  PROBLEM  AND  PROPOSED  SOLUTION  APPROACH 

The  characteristics  of  sensor  clusters  are  different  from  those  of  conventional  ar¬ 
rays  as  given  in  Table  1 .  The  establishment  of  conventional  arrays  is  mostly  localized. 
For  example,  the  antenna  arrays  used  for  radar  and  communications  are  usually  localized 
and  installed  permanently  on  site.  The  beamforming  coordination  is  also  performed  lo¬ 
cally  by  a  central  DSP  processor,  and  the  positions  of  the  elements  and  the  number  of 
elements  are  usually  predetermined  and  fixed  for  a  given  installation.  The  array  elements 
also  remain  operationally  active  for  long  periods  of  time.  On  the  other  hand,  the  estab¬ 
lishment  of  sensor  clusters  is  distributed  in  nature  and  the  positions  of  the  sensor  nodes 
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are  random.  The  number  of  sensor  nodes  required  for  beamforming  is  variable.  The 
situation  is  aggravated  by  the  fact  that  these  sensor  nodes  can  be  mobile  and  subject  to 
frequent  failures  when  their  battery  energy  is  depleted  or  when  they  are  damaged  by  envi¬ 
ronmental  conditions  or  human  intervention. 


No. 

Descriptions 

Conventional  array 

Sensor  cluster  array 

1. 

Array  establishment 

Localized 

Distributed 

2. 

Beamforming  coordination 

Localized 

Distributed 

3. 

Position  of  elements 

Fixed 

Random 

El 

Number  of  elements 

Fixed 

Variable 

5. 

Position  errors  of  elements 

Small 

Can  be  large 

6. 

Element  failures 

Few  or  no  element 

failures 

Frequent  element  fail¬ 
ures 

Table  1 .  Comparison  between  a  conventional  array  and  a  sensor  cluster. 


These  differences  make  the  task  of  distributed  beamforming  using  sensor  clusters 
a  challenging  problem  compared  to  that  of  conventional  arrays,  yet  the  benefits  derived 
from  the  higher  signal-to-noise  ratio  (and  capacity)  gains  due  to  beamforming  are  signifi¬ 
cant  [6]. 

The  proposed  approach  for  distributed  beamforming  using  sensor  cluster  is  de¬ 
scribed  as  follows.  Distributed  array  establishment  and  beamforming  coordination  can  be 
achieved  by  using  clustering  algorithms  to  form  a  hierarchical  clustering  sensor  network 
architecture.  A  primary  node  is  elected  from  among  the  nodes  in  each  cluster  based  on 
some  predefined  criteria.  The  primary  node  is  then  tasked  with  maintaining  frequency, 
phase  and  data  synchronization  among  the  remaining  nodes  (or  secondary  nodes)  within 
the  cluster.  The  primary  node  also  coordinates  the  beamforming  process  to  direct  an  an¬ 
tenna  beam  towards  the  UAV  in  an  adaptive  manner  as  the  UAV  flies  across  the  sensor 
nodes.  The  locations  of  the  sensor  nodes  can  be  determined  using  location  discovery 
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methods,  such  as  trilateration  and  atomic  multilateration  [16].  While  there  are  existing 
methods  to  address  the  distributed  establishment  of  sensor  arrays  and  determination  of 
sensor  node  locations,  there  are  no  existing  methods  to  address  the  number  of  sensor 
nodes  required  for  beamforming  and  the  presence  of  position  errors  and  element  failures. 
Therefore,  in  this  thesis,  the  focus  is  to  tackle  these  issues  by  addressing  the  question  on 
the  number  of  sensor  nodes  required  for  beamforming  and  investigating  the  effects  of  po¬ 
sition  errors  and  element  failures  on  the  quality  of  the  antenna  beam  generated  by  the 
sensor  cluster. 

In  order  to  address  these  issues,  a  simulation  model  was  developed  and  imple¬ 
mented  in  the  MATLAB  programming  language,  and  the  results  of  the  simulation  are 
presented. 


C.  THESIS  OUTLINE 

This  thesis  is  organized  as  follows.  Chapter  II  describes  the  various  WSN  archi¬ 
tectures  and  proposes  a  hierarchical  architecture  for  the  selected  operational  scenario  in 
this  thesis.  This  is  followed  by  the  discussion  of  some  existing  clustering  and  location 
discovery  methods  for  WSN.  The  ideas  of  distributed  beamforming  and  using  sensor 
clusters  as  a  MIMO  system  are  also  described.  Chapter  III  covers  the  fundamental  prin¬ 
ciples  of  antenna  arrays,  random  arrays  and  adaptive  beamforming.  The  use  of  arrays  to 
estimate  the  angle-of-arrival  (AOA)  and  track  the  desired  signal  sent  from  the  UAV  are 
also  discussed.  The  significance  of  the  antenna  beamwidth  on  the  detection  and  tracking 
of  the  UAV  is  highlighted.  Chapter  IV  describes  the  simulation  model  and  presents  the 
simulation  results.  Chapter  V  summarizes  the  work  in  this  thesis  and  offers  several  sug¬ 
gestions  for  future  research.  Finally,  the  Appendix  contains  a  description  of  the 
MATLAB  code  used  in  this  work. 
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II.  INTRODUCTION  TO  WIRELESS  SENSOR  NETWORKS 


This  chapter  discusses  the  various  WSN  architectures  and  proposes  an  architec¬ 
ture  for  the  operational  scenario  used  in  this  thesis.  This  is  followed  by  the  discussion  of 
some  existing  clustering  and  location  discovery  methods.  The  advantages  of  using  sensor 
clusters  to  increase  the  transmission  range  and  the  channel  capacity  between  the  sensor 
nodes  and  the  UAV  are  also  described  in  detail. 

A.  WIRELESS  SENSOR  NETWORK  (WSN)  ARCHITECTURE 

Several  wireless  sensor  network  (WSN)  architectures  have  been  proposed  in  the 
literature.  Figure  3  shows  a  general  WSN  architecture,  which  consists  of  a  sensor  field,  a 
sink  node  and  a  satellite  through  which  the  user  can  access  the  sensor  field  [1].  A  sink 
node  provides  the  link  between  the  satellite  and  the  sensor  nodes,  and  each  sensor  node 
has  the  ability  to  route  the  data  to  the  sink  node.  For  example,  the  data  collected  by  sen¬ 
sor  node  A  is  routed  back  to  the  sink  node  through  node  B ,  C,  D  and  E.  The  sink  node  in 
turn  sends  the  data  to  the  user  node  via  the  satellite. 


Sensor  field 


Figure  3.  WSN  architecture  (After  Ref.  [1].). 


Figure  4  shows  a  two-tiered  hierarchical  clustering  architecture  whereby  sensor 
nodes  are  autonomously  clustered  to  support  energy-efficient  and  scable  operations  [9]. 
A  primary  node  (or  cluster  head)  is  elected  within  each  cluster  to  perform  data  filtering, 
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fusion  and  aggregation.  The  secondary  nodes  communicate  only  with  the  primary  node 
within  each  cluster  while  the  primary  nodes  communicate  with  nodes  in  the  next  tier. 
Sensor  clusters  are  connected  with  each  other  through  gateway  nodes  so  that  data  can  be 
transmitted  to  the  higher  level  nodes  if  no  connectivity  exists  between  a  primary  node  in 
any  cluster  in  Tier  0  and  the  nodes  in  Tier  1 .  A  similar  two-tiered  architecture  has  been 
proposed  in  [10]. 


Tier  1 


C  o 


O- 


...  Sensor  clusters 


TierO 


O  Higher  level  node  @  Gateway  node 
#  Primary  node  Q  Secondary  node 


Figure  4.  Two-tiered  hierarchical  WSN  architecture  (After  Ref.  [9].). 


Figure  5  shows  a  three-tiered  hierarchical  clustering  architecture  proposed  in  [11]. 
It  is  different  from  the  two-tiered  architecture  in  that  there  is  another  tier  on  top  of  Tier  1 
to  consolidate  the  data.  Nodes  in  Tier  1  communicate  with  headquarters  or  a  node  in  Tier 
2.  The  idea  of  a  three-tiered  hierarchical  clustering  architecture  was  also  proposed  by 
[12]  and  [13]. 


Headquarters 


Tier  2 


-O  ... 


O- 


G 


\  Tier  1 


C  o 


o  . 


O  -  •  O 


.-.Sensor  clusters 


•  °  °  )  TierO 


O  Higher  level  node  0  Gateway  node 
0  Primary  node  O  Secondary  node 


Figure  5.  Three-tiered  hierarchical  WSN  architecture  (After  Ref.  [11].). 
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The  various  architectures  presented  are  fairly  similar  to  each  other  in  that  they  are 
hierarchical  architectures.  The  main  advantages  of  tiered  architecture  are  scalability, 
longevity  and  cost-effectiveness  [13].  In  this  thesis,  a  three-tiered  hierarchical  clustering 
WSN  architecture  as  shown  in  Fig.  6  is  considered.  The  sensor  nodes  in  Tier  0  are  being 
combined  to  form  clusters,  and  a  primary  node  is  elected  within  each  cluster.  The  sensor 
clusters  are  connected  with  each  other  through  gateway  nodes  so  that  data  can  be  trans¬ 
mitted  to  the  UAV  if  no  connectivity  exists  between  a  primary  node  in  any  cluster  in  Tier 
0  and  the  UAVs  in  Tier  1 .  The  UAVs  residing  in  Tier  1  are  organized  to  communicate  to 
nodes  at  the  next  higher  tier,  which  are  represented  by  satellite  nodes  in  this  illustration. 
This  aggregation  process  can  be  applied  repeatedly  to  form  a  hierarchy  of  clusters  and  al¬ 
lows  scalable  operation  when  the  number  of  sensor  nodes  or  UAVs  increases  due  to  op¬ 
erational  requirements.  In  this  thesis,  the  focus  is  on  the  communication  link  between  the 
sensor  nodes  and  the  UAV. 


&  Satellite  Tier  2 


@  Gateway  node 
O  Secondary  node 

Figure  6.  Proposed  three-tiered  WSN  architecture. 


The  next  section  describes  some  of  the  clustering  methods  used  to  form  a  hierar¬ 
chical  clustering  WSN  architecture. 
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B.  CLUSTERING  IN  WSN 

The  main  issues  in  the  hierarchical  arrangement  shown  in  Fig.  6  include  how  to 
select  the  primary  node  and  organize  the  cluster.  This  section  discusses  some  of  the 
common  clustering  algorithms  that  researchers  have  proposed  to  address  these  issues. 

1.  Lowest  ID  (LID)  Clustering  Algorithm 

This  is  a  two-hop  clustering  algorithm,  and  Fig.  7  shows  that  a  secondary  node  is 
at  most  two  hops  away  from  another  secondary  node  within  a  cluster.  According  to  [14], 
a  sensor  node  periodically  broadcasts  the  list  of  nodes  (including  itself)  that  it  can  hear. 

A  sensor  node,  which  only  hears  nodes  with  identification  (ID)  numbers  higher  than  itself 
from  its  one-hop  neighborhood,  declares  itself  a  primary  node.  The  newly  declared  pri¬ 
mary  node  then  sets  its  own  ID  number  as  the  cluster  ID  and  broadcasts  this  information 
to  other  nodes.  Nodes  that  are  one-hop  away  from  the  newly  declared  primary  node  re¬ 
ceive  this  information  and  realize  that  they  are  now  secondary  nodes,  leading  to  the  for¬ 
mation  of  a  sensor  cluster. 


Sensor  clusters  are  linked  to  adjacent  clusters  through  assigned  gateway  nodes. 
For  example,  in  Fig.  7,  sensor  node  5  is  a  gateway  node  that  connects  clusters  1  and  3. 
Hence,  a  sensor  node  that  can  hear  two  or  more  primary  nodes  within  two  hops  is  a  gate¬ 
way  node;  otherwise,  it  is  a  secondary  node.  The  main  disadvantage  of  this  method  is 
that  the  energy  dissipation  of  the  sensor  nodes  is  not  evenly  distributed  within  the  sensor 
cluster  since  the  primary  node  is  selected  based  on  ID  numbers,  and  no  other  attributes, 
such  as  node  connectivity  or  energy  level,  are  being  accounted  for. 


Cluster  ID  =  1  Cluster  ID  =  3 


•  Primary  node 
O  Secondary  node 
®  Gateway  node 
Connectivity  exists 


Figure  7.  Lowest  ID  (LID)  clustering  algorithm. 
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2.  Highest  Connectivity  (HCN)  Clustering  Algorithm 

The  highest  connectivity  (HCN)  clustering  algorithm  is  similar  to  the  LID  algo¬ 
rithm  except  that  a  node  is  elected  as  the  primary  node  if  it  has  the  largest  number  of 
connected  nodes  around  it  [14].  As  shown  in  Fig.  8,  nodes  1  and  3  are  elected  as  primary 
nodes  because  both  of  them  have  the  largest  number  of  connected  nodes.  This  algorithm 
is  better  than  the  LID  algorithm  in  that  HCN  takes  connectivity  into  account. 

•  Primary  node 
O  Secondary  node 
@  Gateway  node 
—  Connectivity  exists 


Figure  8.  Highest  connectivity  (HCN)  clustering  algorithm. 


3.  Low-energy  Adaptive  Clustering  Hierarchy  (LEACH) 

This  algorithm  aims  to  minimize  the  energy  dissipation  in  WSN  by  changing  the 
primary  node  assignment  periodically  so  that  sensor  nodes  take  turns  to  perform  the  role 
of  a  primary  node  for  a  certain  period  of  time  defined  as  a  round.  Since  data  transfer  is 
an  energy  consuming  process  and  is  done  by  the  primary  node,  this  rotation  of  primary 
nodes  leads  to  more  balanced  energy  dissipation  in  the  entire  network  compared  to  the 
LID  and  HCN  algorithms. 

The  LEACH  algorithm  consists  of  two  phases  [15]:  setup  and  steady.  During  the 
setup  phase,  each  sensor  node  chooses  a  uniform  random  number  between  0  and  1  and  if 
the  number  is  lower  than  the  threshold  T(n)  for  node  /?,  the  sensor  node  becomes  the  pri¬ 
mary  node.  The  threshold  T(n)  for  node  n  is  given  by  [15] 


P 


T{n)  =  { 


l-P[/*mod(l/P)] 


0 


if  n  e  G 
otherwise. 


(i) 


ii 


where  P  is  the  desired  percentage  of  sensor  nodes  that  are  expected  to  be  primary  nodes, 
r  is  the  current  round,  G  is  the  set  of  sensor  nodes  that  has  not  been  primary  nodes  in  the 
past  1/P  rounds  and  mod( )  denotes  the  modulus  after  division.  Equation  (1)  is  used  to 
ensure  that  all  sensor  nodes  dissipate  equal  amount  of  energy.  Once  the  primary  node  has 
been  selected,  the  next  phase  is  the  steady  phase  in  which  the  primary  node  assigns  a 
time-division  multiple  access  (TDMA)  schedule  to  the  secondary  nodes  so  that  the  latter 
can  transfer  their  data  to  the  primary  node  at  predefined  time  slots.  The  steady  phase 
lasts  for  a  longer  duration  than  the  setup  phase.  After  a  certain  period  of  time  in  the 
steady  phase,  the  setup  phase  is  initiated  again,  leading  to  a  change  in  the  primary  node 
assignment. 

In  summary,  three  types  of  clustering  algorithms  were  introduced  in  this  section: 
primary  node  selection  based  on  ID  numbers,  highest  connectivity  and  energy  dissipa¬ 
tion.  Now,  we  proceed  to  discuss  the  task  of  determining  the  location  of  the  sensor  nodes 
in  the  next  section. 

C.  LOCATION  DISCOVERY  IN  WSN 

The  objective  of  deploying  a  WSN  is  to  detect  the  presence  of  targets  of  interest 
and  hence  identify  their  positions.  This,  in  turn,  poses  a  requirement  on  the  sensor  nodes 
to  determine  their  own  positions,  thus  the  need  for  location  discovery.  This  section  dis¬ 
cusses  some  of  the  location  discovery  techniques  that  can  be  used  to  determine  the  loca¬ 
tions  of  the  sensor  nodes,  assuming  that  one  or  more  beacon  nodes  exist  [16].  Beacon 
nodes  are  nodes  that  have  absolute  information  of  their  own  location.  Generally,  this  as¬ 
sumption  is  valid  because  in  WSN,  only  a  small  number  of  sensor  nodes  may  be 
equipped  with  GPS  receivers  due  to  economic  reasons.  They  function  as  location  refer¬ 
ences  by  periodically  broadcasting  their  location  and  hence  allowing  other  sensor  nodes 
(not  equipped  with  GPS  receivers)  to  determine  their  locations. 

The  common  methods  of  measuring  relative  distances  between  beacon  nodes  and 
a  given  node  are  listed  as  followed  [16]: 

1 .  A  given  sensor  node  A  may  use  the  time  of  arrival  (ToA)  techniques  to 
measure  the  propagation  time  of  a  signal  emitted  by  a  beacon  node  B  and 
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then  calculate  the  relative  distance  between  itself  and  the  beacon  node  as¬ 
suming  that  the  speed  of  the  signal  in  the  medium  is  known. 

2.  Time  difference  of  arrival  (TDoA)  techniques  utilize  two  or  more  simulta¬ 
neously  transmitted  signals  from  different  beacon  nodes  to  determine  the 
time  difference  between  them.  The  distance  from  the  sensor  node  to  a  bea¬ 
con  node  can  be  accurately  calculated  using  the  measured  TDoA. 

3.  A  sensor  node  may  also  utilize  the  received  signal  strength  indicator  (RSSI) 
to  determine  the  strength  of  a  received  signal  and  hence  calculate  the  dis¬ 
tance  between  itself  and  the  beacon  node  assuming  that  the  propagation  loss 
between  them  is  known. 


After  the  distance  between  the  sensor  node  and  the  beacon  nodes  are  determined 
based  on  the  methods  described  above,  there  remains  a  need  to  combine  these  distance 
measurements  to  find  the  actual  location  of  the  sensor  node.  According  to  [16],  the  trilat- 
eration  and  atomic  multilateration  methods  are  the  two  most  common  ways  to  combine 
the  distance  measurements  from  three  or  more  beacons.  The  trilateration  and  atomic 
multilateration  methods  are  illustrated  in  Fig.  9  and  described  in  the  following. 


Trilateration 


Atomic  multilateration 


Beacon  node  with  known  location 
Sensor  node  with  unkown  location 
Distance  measurement 


Figure  9.  Trilateration  and  atomic  multilateration  for  node  location  discovery 

(After  Ref.  [16].). 


Trilateration  is  used  when  a  node  has  accurate  distance  measurements  between  it¬ 
self  and  at  least  three  beacon  nodes  [16].  In  Fig.  9,  the  intersection  of  the  three  circles 
each  centered  at  beacons  Bj ,  and  Bs  gives  the  location  of  the  sensor  node  A.  How- 
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ever,  this  method  may  not  work  well  if  the  circles  do  not  intersect  at  a  single  point  [16]. 
On  the  other  hand,  atomic  multilateration  uses  more  than  three  beacons  to  discover  the 
location  of  a  sensor  node.  According  to  [16],  this  technique  defines  an  objective  function 
that  minimizes  the  error  between  the  different  measured  distances  by  calculating  a  local 
minimum  of  that  function.  Consider  the  case  in  which  the  exact  location  of  node  A  is  de¬ 
noted  by  A(X,Y),  and  the  locations  of  beacon  nodes  are  denoted  by  B^(Xx,Yx ),..., 

BK  ( XK ,  Yk  ).  One  type  of  objective  function  is  the  sum  of  residuals  as  defined  by  [  1 6] 

Ll(X,Y)  =  t\Di(X,Y)-Ri(X,Y)\  (2) 

i=l 

where  the  distance  D.(X,  Y)  is  given  by 

Di(X,Y)  =  J{X,-X?  +  {Yi-Yj‘  (3) 

and  R.(X,Y)  is  the  estimated  distance  between  beacon  node  B. (X. , Y. )  and  the  node 
A(X ,  Y).  A  local  minimum  of  the  function  given  by  Equation  (2)  is  computed  in  order  to 
determine  the  location  of  node  A. 

D.  TRANSMISSION  RANGE  EXTENSION 

After  the  clustering  and  location  discovery  functions  are  accomplished,  there  re¬ 
mains  the  question  of  transmitting  the  data  back  to  the  UAV.  As  shown  in  Fig.  6,  the  role 
of  data  transmission  is  taken  up  by  the  primary  node  in  each  sensor  cluster.  However,  the 
limited  energy  available  within  the  primary  node  and  the  omni-directional  nature  of  its 
antenna  radiation  pattern  may  restrict  the  transmission  range  between  itself  and  the  UAV. 
Therefore,  it  is  proposed  in  this  thesis  to  use  a  sensor  cluster  to  form  an  array-like  trans¬ 
mission  to  the  UAV. 

Figure  10  illustrates  the  potential  additional  antenna  gain  G\  that  can  be  achieved 
by  using  a  sensor  array  compared  to  that  of  a  single  sensor  node.  An  additional  gain  Gj 
may  be  achieved  by  using  more  sensor  nodes. 
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Antenna  beam  generated  Antenna  beam  generated  by  Antenna  beam  generated  by 
by  a  sensor  node  a  sensor  cluster  with  N-  nodes  a  sensor  cluster  with  Ns »  N- 


Figure  10. 


Higher  antenna  gain  derived  from  beamforming  using  sensor  cluster. 


The  free  space  transmission  loss  can  be  expressed  as  [3] 

P  _  G.G, 

P  K 

where  P  is  the  power  received  by  the  receiver,  Pt  is  the  power  transmitted  by  the  trans¬ 
mitter,  G  is  the  antenna  gain  of  the  receiver,  Gt  is  the  antenna  gain  of  the  transmitter, 
and  Ls  is  the  free  space  loss  given  by 

L.-^f  (5) 

where  R  is  the  transmission  range  and  A,  is  the  free  space  wavelength.  If  the  sensor 
cluster  and  the  UAV  are  represented  as  transmitter  and  receiver,  respectively,  it  can  be 
shown  in  Equation  (4)  that  for  fixed  Pr/P,  ratio,  an  increase  in  G,  by  a  factor  of  y  can 
lead  to  a  reduction  in  L  by  the  same  factor.  This  in  turn  increases  the  transmission  range 
R  by  a  factor  of  yjy  using  Equation  (5).  Hence,  a  sensor  cluster  can  be  used  to  form  an 
antenna  beam  with  higher  gain  to  increase  the  transmission  range  for  data  transmission. 
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E.  CHANNEL  CAPACITY  EXPANSION 

The  data  rate  capacity  between  the  sensor  nodes  and  the  UAV  can  be  increased  by 
exploiting  the  multiple-input  multiple-output  (MIMO)  communication  nature  of  the  sys¬ 
tem,  which  can  provide  data  rate  gains  [4].  MIMO  systems  attempt  to  create  multiple 
transmission  channels  using  several  independent  mulipaths  in  a  multipath  rich  environ¬ 
ment  to  increase  the  channel  capacities  as  compared  to  a  single  channel  using  a  transmit¬ 
ting  and  receiving  antenna  pair.  A  basic  two-input-two-output  narrowband  MIMO  sys¬ 
tem  is  shown  in  Fig.  11  where  four  transmission  paths  exist  between  the  transmitter  and 
the  receiver  antenna  arrays. 

T  ransmitter 

P  i,  n  ,  PH,  nH 

Figure  1 1.  Two-input-two-output  narrowband  MIMO  system. 

Let  nr  and  nR  be  the  number  of  antenna  elements  at  the  transmitter  and  receiver, 
respectively.  The  total  transmitted  power  is  indicated  by  Pn  and  PR  is  the  average  power 
over  all  receive  elements.  The  average  signal-to-noise  ratio  at  each  receiving  antenna 
branch  is  then  given  by  [4] 

P  =  4  (6) 

CJ~ 

where  a2  is  the  noise  power  at  each  receiving  antenna  branch. 

The  channel  impulse  response  matrix  G  is  given  by 


where  it  is  assumed  that  the  channel  impulse  response  is  constant  over  the  frequency 
band  of  interest,  hence  the  term  “narrowband  MIMO  system”.  The  commonly  used  nor¬ 
malized  channel  impulse  response  matrix  is  given  by 
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H  = 


'  P\ 
4L  G 

,PJ- 


(8) 


and  the  capacity  of  the  overall  MIMO  channel  C  can  be  represented  as  the  sum  of  the 
classical  Shannon  capacities  of  each  of  the  spatial  channels  given  by  [4] 


C  =  log2  det 


"R 


nT 


(9) 


where  det  [•]  represents  the  determinant  and  "  represents  the  conjugate  transposition.  Let 
Xk  represent  the  eigenvalues  of  HH" .  In  the  case  of  the  two-input-two-output  narrow- 
band  MIMO  system,  Equation  (9)  is  equivalent  to  the  following  [4] 


C  =  Z  log2 

A=1 


(10) 


which  is  simplified  to  become 


C  =  Iog2 


i+f^ 


log 


1  +— A, 


(11) 


Assuming  that  the  channel  is  lossless  and  all  eigenvalues  are  the  same,  it  can  be 
shown  in  Equation  (11)  that  the  channel  capacity  increases  by  a  factor  of  two.  Figure  12 
shows  a  general  /?7  -input-/?,,  -output  narrowband  MIMO  system  with  the  normalized 
channel  impulse  response  given  by  Hir  Considering  the  increased  number  of  transmis¬ 
sion  multipaths  between  the  transmitter  and  receiver  as  shown  in  Fig.  12,  Equation  (10) 
can  be  extended  so  that  the  capacity  C  becomes 


C  =  Z  log, 

A=l 


(12) 


where  K  is  mm[nT,nR].  Therefore,  the  channel  capacity  can  be  increased  to  several  or¬ 
ders  of  magnitude  by  increasing  n,  and  nR,  as  indicated  by  Equation  (12). 
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Pr,nT 


Pr>  Hr 


Figure  12.  A  MIMO  narrowband  system  with  nT  inputs  and  nt  outputs. 


The  sensor  nodes  may  be  deployed  over  an  area  where  a  large  number  of  commu¬ 
nication  transmission  multipaths  exist  due  to  the  presence  of  man-made  structures  or 
natural  obstacles.  Figure  13  illustrates  how  the  sensor  cluster  and  the  UAV  can  be  con¬ 
sidered  a  type  of  MIMO  system  consisting  of  multi-element  array  elements  at  both  ends 
of  the  communication  link.  Multiple  channels  are  created  using  several  independent  mul¬ 
tipaths  in  the  multipath  rich  environment  to  increase  the  channel  capacities.  Such  a  sys¬ 
tem  could  offer  a  unique  solution  for  supporting  large  bandwidth  requirements  in  future 
military  applications. 


Sensor  cluster 


Figure  13.  Sensor  cluster  as  part  of  MIMO  system. 
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F.  SUMMARY 

In  this  chapter,  various  WSN  architectures  were  discussed,  and  a  hierarchical 
three-tier  architecture  is  proposed  for  use  in  this  thesis.  This  was  followed  by  the  discus¬ 
sion  of  some  common  clustering  and  location  discovery  methods  for  WSNs.  The  advan¬ 
tages  of  using  sensor  cluster  to  increase  the  transmission  range  and/or  channel  capacity 
between  the  sensor  nodes  and  the  UAV  were  also  described.  The  next  chapter  will  dis¬ 
cuss  the  fundamental  concepts  required  for  each  sensor  cluster  to  form  an  array-like 
transmission  to  the  UAV. 


19 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


20 


III.  SENSOR  CLUSTER  AND  BEAMFORMING 


After  the  sensor  nodes  are  deployed  by  the  UAV,  they  are  grouped  into  clusters, 
and  a  primary  node  is  elected  from  among  the  sensor  nodes  within  a  cluster  using  the 
methods  described  in  Chapter  II,  Section  B.  The  location  of  the  sensor  nodes  can  be  de¬ 
termined  using  techniques  presented  in  Chapter  II,  Section  C.  Distributed  beamforming 
is  then  initiated  and  coordinated  by  each  primary  node  within  the  cluster  in  order  to  in¬ 
crease  the  transmission  range  and/or  channel  capacity  between  each  cluster  and  the  UAV. 

This  chapter  presents  the  fundamental  concepts  required  for  each  sensor  cluster  to 
form  an  array-like  transmission  to  the  UAV.  In  particular,  it  discusses  the  array  factor 
and  power  gain  of  antenna  arrays  and  random  arrays.  The  adaptive  beamforming  scheme 
is  also  described.  This  is  followed  by  the  discussion  of  various  algorithms  that  can  esti¬ 
mate  the  AOA  and  track  the  desired  signal  using  passive  arrays.  Lastly,  it  also  discusses 
the  significance  of  the  antenna  beamwidth  for  UAV  detection  and  tracking. 


A.  ONE-DIMENSIONAL  ARRAY 

Figure  14  shows  an  A/xl  omni-directional  antenna  array  located  in  the  far  field  of 
a  point  source  (e.g.,  an  UAV).  The  point  source  emits  a  signal  that  consists  of  a  complex 
signal  s(t)  modulating  a  carrier  e"”'\  The  plane  wavefront  arrives  at  an  elevation  angle 
0{)  with  respect  to  the  array  normal  (i.e.,  z-axis).  For  convenience,  the  reference  is  taken 


at  the  origin,  and  it  is  assumed  that  the  wavefront  arrives  at  the  m'h  element  tm  seconds  be¬ 
fore  it  reaches  the  origin.  Hence,  the  signal  arriving  at  the  m'h  element  leads  the  signal  ar¬ 
riving  at  the  origin  by  [5] 


UA) 


c 


(13) 


where  c  is  the  speed  of  light  and  xm  is  the  A'-coordinate  of  the  mth  element. 
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Array  Normal 
z-axis 


Figure  14.  An  Mxl  omni-directional  antenna  array. 


It  can  be  seen  from  Equation  (13)  that  tm(0Q )  is  determined  by  the  AOA  of  the  de¬ 
sired  signal  indicated  by  and  the  element’s  position  indicated  by  xm.  For  any  arbitrary 
angle  d ,  the  signal  arriving  at  the  m'h  element  leads  the  signal  arriving  at  the  origin  by 


Without  considering  the  effects  of  mutual  coupling  among  the  array  elements,  the 
spatial  response  (or  array  factor)  of  the  array  is  obtained  by  adding  all  the  array  element 
outputs  together  giving  [  1 8] 


F{6)  =  I  Ie 


-iOc'mWgi'Oc'mW) 


where  Im  and  e~,<0c'mW  are  the  magnitude  and  phase  of  the  current  induced  on  the  m'h  ele¬ 
ment.  Together,  Im  and  e~'<0c,M))  form  what  is  commonly  known  as  the  complex 
weights.  Notice  that  the  phase  reversal  in  the  term  is  required  in  order  to  create 


a  maximum  value  of  F(0)  at  0  =  0{).  Hence,  the  maximum  value  of  F(0)  occurs  at 
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9  =  9Q  and  the  main  lobe  points  towards  0Q.  Equation  (15)  can  be  written  in  the  form 
given  by  [18] 


F(0)  =  ll„e 

m= 1 


sin%  >e  )p{xm  sin  0) 


(16) 


where  fi  is  the  wavenumber  given  by  2 n!  X  and  X  is  the  wavelength,  and  can  be  com¬ 
pactly  expressed  as 


F(0)  =  £  we 


jp(xmsmO) 


(17) 


where  w,„=hemxM  is  the  complex  weight  applied  to  the  ni"'  element.  The  far-field 
beam  pattern  (or  power  gain)  is  given  by 

G{0)  =1  F(0)  |2 .  (18) 

Consider  an  A/xl  uniformly  excited  unit  array  (i.e.,  equal  current  magnitude  of  1) 
consisting  of  isotropic  elements  (i.e.,  antennas  with  omni-directional  radiation  pattern) 
and  with  the  main  lobe  pointing  at  90  =  0  .  The  maximum  array  factor  is  computed  to  be 


m)=s>> 


'*)  - 


=  1.1  e 


-jP(xmsxn0Q) ^jp{xmsm0O)  _  ^  j  _ 


using  Equation  (17)  while  the  power  gain  is  calculated  to  be  M~  using  Equation  (18). 
Therefore,  the  power  gain  of  a  uniform  array  can  be  increased  by  increasing  the  number 
of  elements.  This  implies  that  more  sensor  nodes  can  be  used  within  a  sensor  cluster  to 
produce  an  antenna  beam  with  higher  power  gain  and  hence  increase  the  transmission 
range  between  the  sensor  cluster  and  the  UAV,  which  was  described  earlier  in  Chapter  II, 
Section  D.l. 

The  plot  of  the  normalized  power  gain  Gn  (in  dB)  versus  9  (in  degrees)  is  shown 
in  Fig.  15  to  illustrate  the  main  characteristics  of  the  antenna  beam  generated  by  a  10x1 
uniformly  excited  array.  The  array  consists  of  isotropic  elements  with  fixed  element- 
element  spacing  d  =  A/2,  and  the  main  lobe  is  assumed  to  point  at  90  =  0  Figure  15(a) 

shows  that  the  main  lobe  points  in  a  direction  perpendicular  to  the  array  axis  (i.e.,  x-axis) 
and  can  be  categorized  as  a  broadside  radiation  pattern.  As  shown  in  Fig.  15(b),  the  3-dB 


23 


0 
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Normalized  power  gain,  G„  (dll) 


(a) 


0  (degrees) 


(b) 

Figure  15.  Normalized  power  gain  of  the  beam  generated  by  a  10x1  uniformly  excited 
square  array  with  isotropic  elements,  fixed  element  spacing  of  d  =  X/2,  and  0O  =  0  :  (a) 

Polar  plot  and  (b)  X-Y  plot  indicating  the  3-dB  beamwidth  of  10.2  and  the  highest 

sidelobe  level  of  -13.2  dB. 
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(or  half  power)  beamwidth  #3dB  is  about  10.2°,  and  the  highest  sidelobe  level  SLmax  is 
about  13.2  dB below  the  main  lobe  level. 


The  3-dB  beamwidth  of  a  uniformly  excited  broadside  array  (i.e.,  6{)  =  0  )  with 
isotropic  elements  and  fixed  element  spacing  is  given  by  [18] 


ft  ««  0.866- 


X 


3dB 


Md 


(19) 


where  a  long  (i.e.,  Md  »  X)  array  has  been  assumed.  The  3-dB  beamwidth  is  calculated 

to  be  9.92  using  Equation  (19),  showing  reasonable  agreement  with  the  10.2  obtained 
from  the  plot.  It  can  be  seen  from  Equation  (19)  that  ftdB  can  be  controlled  by  adjusting 
M.  For  example,  when  M  is  increased  by  a  factor  of  y,  #?dB  is  reduced  by  the  same  fac¬ 
tor,  resulting  in  a  narrower  beam  and  vice  versa.  This  allows  03dB  to  be  adjusted  by 
changing  the  number  of  sensor  nodes  required  for  beamforming  during  the  detection  and 
tracking  of  UAV.  For  example,  a  large  #3dB  is  usually  required  to  create  a  large  detection 
region  while  a  small  #3dB  is  preferred  to  obtain  the  required  tracking  resolution  of  the 
UAV.  This  is  further  discussed  in  Section  G.2  later  in  this  chapter. 

The  array  factor  and  the  power  gain  of  a  one-dimensional  array  can  be  extended 
to  a  two-dimensional  array  as  presented  in  the  next  section. 


B.  TWO-DIMNESIONAL  ARRAY 

Figure  16  shows  an  MxN  omni-directional  antenna  array  located  in  the  far  field  of 
a  point  source.  The  plane  wavefront  arrives  at  an  elevation  angle  ft0  with  respect  to  the  z- 
axis  and  an  azimuth  angle  of  j>0  with  respect  to  the  x-axis.  The  reference  is  taken  at  the 
origin,  and  it  is  assumed  that  the  wavefront  arrives  at  the  (m,n)'h  element  tmn  seconds  be¬ 
fore  it  reaches  the  origin.  Therefore,  the  signal  arriving  at  the  ( m,r)'h  element  leads  the 
signal  arriving  at  the  origin  by  [5] 

,  «  .  .  .x;  sin  ft  cos +  y  sin  ft  sin  <An 

L,  (e„  A,)= — — - — — - — -  (20) 

C 
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where  c  is  the  speed  of  light  and  (xmn,ymn)  is  the  x-coordinate  and  ^-coordinate  of  the 
( m,n)lh  element. 


Array  Normal 
z-axis 


Array  Normal 
z-axis 


Point  source,  s(t)e,<ac' 


Figure  16.  An  MxN  omni-directional  antenna  array. 


For  any  set  of  arbitrary  angles,  (6,j)) ,  the  signal  arriving  at  the  (m,/?)lh  element 
leads  the  signal  arriving  at  the  origin  by 

xmn  sin  6  cos  (f>  +  ymn  sin  6 sin  ^ 


tA0,4)  = 


(21) 


The  two-dimensional  array  factor  F(0,(f)  is  given  by  [18] 


M  N 


mi)  =11.1  e 


-  iPixmn  sin 4,  cos ^  +ymn  sin  sin  ^ )  jfi( xmn  sin 0cosf+ymn  sin  0  sin**) 


(22) 


m= I  «=1 


without  considering  the  effects  of  mutual  coupling  among  the  array  elements.  Equation 
(22)  can  be  rewritten  as 


M  N 


m,i)= isv 


sin^cos^+ymn  sin  £  sin**) 


(23) 


where  /?  is  the  wavenumber  given  by  2 n!  k,  and  wmn  =  I mnejP{Xmn  s,n ^ cos ^ s,n ^ s,n ^ >  ;s  the 
complex  weight  applied  to  the  (m,n)"'  element.  The  maximum  value  of  F(6,j))  occurs 
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at  (#,^)  =  (0O,^O),  and  the  main  lobe  points  towards  (0O,^O).  The  power  gain  is  defined 
as 

G(0,^)=\F(0J)\2 .  (24) 

Consider  an  MxN  uniformly  excited  array  with  isotropic  elements  and  the  main 
lobe  pointing  at  0Q  =  </0  =  0 .  The  maximum  array  factor  and  the  power  gain  are  NM  and 

[NM)~ ,  respectively.  Similar  to  the  one-dimensional  case,  the  power  gain  of  a  two  di¬ 
mensional  uniform  array  can  be  increased  by  increasing  the  number  of  elements. 

The  beam  pattern  formed  by  3x3,  5x5  and  7x7  arrays  are  shown  in  Fig.  17  in 
which  plots  were  generated  for  (f  =  45  .  These  uniformly  excited  square  arrays  consist  of 
isotropic  elements  with  a  fixed  element-element  spacing  d  =  )J2.  The  main  lobe  for  all 
the  three  arrays  is  assumed  to  point  at  (#0,^0)  =  (0  ,45  )  as  illustrated  in  Fig.  17(a). 
Figure  17(b)  shows  that  #JdB  increases  from  15  to  37  when  the  number  of  array  ele¬ 
ments  decreases  from  7x7  to  3x3.  Therefore,  the  antenna  beam  of  a  uniformly  excited 
array  with  isotropic  elements  and  fixed  element-element  spacing  becomes  broader  when 
the  number  of  elements  decreases.  This  observation  is  consistent  with  that  concluded  in 
the  previous  section  for  the  one-dimensional  case. 
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7x7  array  - 5x5  array  - 3x3  array  2-D  plots  generated  at  4*  =  45' 


Figure  17.  Normalized  power  gain  of  the  beams  generated  by  7x7,  5x5  and  3x3  uniformly 
excited  square  arrays  with  isotropic  elements,  fixed  element  spacing  of  d  =  A./2,  and 
(#0,^0)  =  (0  ,45  ) :  (a)  Polar  plots  and  (b)  X-Y  plots  showing  the  decrease  in  3-dB  beam- 

width  as  the  number  of  elements  increase. 
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In  summary,  the  array  factor  and  power  gain  equations  for  the  two-dimensional 
array  were  presented  in  this  section,  and  it  was  concluded  that  they  were  affected  by  the 
AOA  of  the  desired  signal,  the  number  of  array  elements,  array  element  positions,  and  the 
complex  weights  wmn.  The  main  characteristics  of  the  antenna  beam,  such  as  3-dB 

beamwidth  and  highest  sidelobe  level,  were  introduced.  The  effects  of  the  number  of  ar¬ 
ray  elements  on  the  antenna  power  gain  and  the  3-dB  beamwidth  were  also  discussed. 

C.  RANDOM  ARRAY 

Sensor  nodes  are  deployed  randomly  over  an  area  of  interest,  resulting  in  random 
positioning  of  the  sensor  nodes.  Subsequently,  these  sensor  nodes  form  clusters  using 
methods  described  in  Chapter  II,  Section  B  and  their  locations  are  determined  using  loca¬ 
tion  discovery  techniques  presented  in  Chapter  II,  Section  C.  The  random  position  of  the 
sensor  nodes  within  the  sensor  cluster  gives  rise  to  what  is  referred  to  as  random  arrays. 
This  section  discusses  the  theoretical  findings  related  to  random  arrays,  positioning  error 
and  element  failures  in  antenna  arrays. 


1.  Average  Power  Gain 

The  average  power  gain  of  one  dimensional  A/xl  random  array  is  given  by  [7] 


G(0)  =  \F(9) 


1  (1| 

M  )\M  j 


(25) 


where  Go  (#)  is  the  designed  power  gain  and  M  is  the  number  of  array  elements. 


The  first  term  Go  (0)^1  -  —J  is  an  angle-dependent  term,  which  shows  that 

Go  ( 6 )  is  reduced  depending  on  the  number  of  array  elements  M.  For  example  M  must 
decrease  from  100  to  2  in  order  to  produce  a  3-dB  reduction  in  G0  (0).  The  second  term 

is  an  angle-independent,  additive  term  and  relates  to  the  power  ratio  of  the  average 

sidelobe  to  the  main  lobe.  For  example,  the  average  sidelobe  level  is  computed  to  be 
-20  dB  when  M  =  100  and  -3  dB  when  M  =  2,  an  increment  of  1 7  dB.  It  can  be  seen 
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that  increasing  the  number  of  array  elements  can  help  reduce  the  loss  in  the  average 
power  gain  and  decrease  the  average  sidelobe  level,  even  though  the  latter  is  more  sensi¬ 
tive  to  the  number  of  elements  than  the  average  power  gain.  It  was  further  suggested  in 
[7]  that  the  shape  of  the  main  lobe  is  also  not  significantly  affected  by  the  random  posi¬ 
tions  of  the  array  elements. 

2.  Position  Errors 

Reference  [8]  considered  the  presence  of  both  amplitude  and  phase  errors  in  a 
two-dimensional  MxN  array  with  equal  weights,  i.e.,  =  /.  The  expected  value  of  the 

increase  in  sidelobe  level  As,  relative  to  the  main  lobe,  is  given  by 


A.9  = 


MN 


1  +  %- 

i-  j 


(26) 


where  cr\a  is  the  variance  of  the  amplitude  error  in  volts2  and  <j 2N,  is  the  variance  of  the 
phase  error  in  rad2,  assuming  gaussian  phase  errors. 

Considering  only  the  presence  of  phase  errors,  it  can  be  calculated,  using  Equa¬ 
tion  (26),  that  the  sidelobe  level  increases  about  6  dB  when  cr^  increases  by  a  factor  of 

two.  For  example,  assuming  that  MN  =  9,  As  is  calculated  to  be  0.0034  when  <j^=  10 
and  As  =  0.014  when  a 20",  an  increment  of  about  6  dB  caused  by  an  increase  in  the 
phase  error  by  a  factor  of  two. 

It  is  also  noted  from  Equation  (26)  that  the  amount  of  increase  in  the  sidelobe 
level  can  be  reduced  by  using  more  array  elements  indicated  by  the  MN  term  in  the  de¬ 
nominator.  It  can  be  shown  that  it  takes  the  number  of  array  elements  to  increase  by  a 
factor  of  two  in  order  to  reduce  As  by  3  dB. 

The  fractional  loss  Lp  in  the  main  lobe  gain  for  a  two-dimensional  MxN  array  with 
equal  weights  is  given  by  [8] 

L,  =  e°%’  (27) 


where  is  the  variance  of  the  phase  error,  assuming  gaussian  phase  errors.  It  can  be 

calculated,  using  Equation  (27),  that  a  phase  error  of  ±47  is  required  to  cause  a  3-dB  re- 
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duction  in  the  main  lobe  gain.  This  suggests  that  a  large  position  error  is  required  to  re¬ 
duce  the  main  lobe  gain  by  3  dB. 

3.  Element  Failures 

The  fractional  loss  L/  in  the  main  lobe  gain  for  a  two-dimensional  MxN  array  with 
equal  weights  due  to  element  failures  is  given  by  [8] 


mn'  y 

MN  ) 


(28) 


where  MW' is  the  number  of  operating  elements  and  MN  is  the  total  number  of  ele¬ 
ments.  Since  element  failures  are  random,  the  operating  probability  of  any  element  p0  is 
given  by 


Po 


MN' 
Inn - 

A/V— *ao  j^N 


(29) 


Hence,  the  expected  value  of  the  main  lobe  gain,  relative  to  error-free  array,  is  re¬ 
duced  by  a  factor  of  pi . 


In  summary,  this  section  has  developed  some  insights  into  the  effects  of  random 
positions  of  array  elements,  random  position  errors  and  element  failures  on  the  antenna 
beam  characteristics.  Some  important  conclusions  can  be  made  as  follows: 

1 .  The  average  power  gain  of  the  main  lobe  is  not  significantly  affected  by  random  posi¬ 
tions  or  position  errors  of  the  array  elements.  The  shape  of  the  main  lobe  is  also  not 
significantly  affected  by  the  random  positions  of  the  array  elements. 

2.  Sidelobe  levels  are  affected  by  random  positions  or  position  errors  of  the  array  ele¬ 
ments,  and  a  lower  sidelobe  level  is  likely  to  be  maintained  when  the  number  of  array 
elements  is  large. 

3.  The  loss  in  the  main  lobe  power  gain  is  proportional  to  the  square  of  the  probability 
of  the  element  being  operative. 

After  the  fundamental  ideas  of  two  dimensional  arrays  and  random  arrays  were 
presented,  the  next  section  describes  the  array  beamforming  processing  which  is  required 
to  form  a  directed  antenna  beam  towards  the  UAV. 
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D.  BEAMFORMING 

It  was  concluded  in  Section  B  of  this  chapter  that  the  array  factor  is  affected  by 
the  AOA  of  the  desired  signal,  the  number  of  array  elements,  array  element  positions,  and 
the  complex  weights  wmn.  An  array  beamformer  works  by  adjusting  the  magnitudes 
and/or  phases  of  the  complex  weights  wnw  and  steers  the  main  lobe  towards  the  desired  di¬ 
rection,  assuming  that  the  AOA  of  desired  signal,  the  number  of  array  elements  and  posi¬ 
tions  of  the  array  elements  are  known.  In  a  WSN,  the  positions  of  the  sensor  nodes  can 
be  obtained  by  methods  described  in  Chapter  II,  Section  C  while  methods  to  estimate  and 
track  the  AOA  of  the  desired  signal  are  described  in  Section  F  of  this  chapter.  This  sec¬ 
tion  starts  by  discussing  some  of  the  common  array  beamformers  found  in  [17]  before 
presenting  the  type  of  beamformer  selected  for  this  thesis. 

According  to  [17],  a  conventional  beamformer  is  a  delay-and-sum  beamformer 
with  complex  weights  of  equal  magnitudes.  The  phases  are  selected  so  as  to  steer  the 
beam  in  the  desired  direction.  Another  type  of  beamformer  called  the  null-steering  beam- 
former  adjusts  the  complex  weights  so  as  to  produce  a  null  in  the  array  factor  in  the  direc¬ 
tion  of  the  interference  signal.  An  optimal  beamformer  computes  the  complex  weights  so 
that  the  main  lobe  has  the  highest  signal-to-noise  ratio.  A  narrowband  beamformer  is  ba¬ 
sically  a  combination  of  the  conventional  and  null-steering  beamformers.  It  uses  a 
known  reference  signal  (which  is  closely  correlated  to  the  desired  signal)  to  estimate  the 
complex  weights  by  subtracting  the  array  output  from  this  reference  signal  and  comput¬ 
ing  the  error.  The  complex  weights  are  then  adjusted  so  as  to  minimize  the  mean  square 
error.  This  type  of  beamformer  has  the  ability  to  steer  the  beam  in  the  desired  direction 
and  create  nulls  in  the  direction  of  the  interfering  signals.  As  such,  the  narrowband 
beamformer  is  used  in  this  thesis. 

Consider  an  MxN  omni-directional  antenna  array  shown  in  Fig.  1 8.  Recall  from 
Fig.  16  that  the  point  source  emits  a  signal  that  consists  of  a  complex  modulating  signal 
s(t)  modulated  by  a  carrier  eio>c' .  The  signal  received  at  the  (/;?,/?)Ih  element  at  time  t  is 
given  by 

*  „  =  +  <j>0 ))  +  nmn  (/)  (30) 
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where  tmn  (&09</>0)  represents  the  signal  time  lead  at  the  (m,n)'b  element  with  respect  to  the 
origin  and  nmn(t)  is  white  noise  with  zero  mean  and  variance  of  cr.  As  shown  in  Fig.  18, 
the  complex  weights  wmn,  which  were  introduced  in  Section  B,  are  applied  to  the  signal 
received  by  the  array  elements  and  summed  together  to  give  the  beamformer  output. 


Figure  18.  MxN  narrow-band  beamformer. 

Generally,  s(t)  is  assumed  to  be  a  narrow-band  signal  and  hence  satisfies  the  fol¬ 
lowing  condition  given  by  [19] 

BAT  «1  (31) 

where  Bm  is  the  bandwidth  of  s(t)  and  Ariax  is  the  maximum  travel  time  of  s(t)  between 
two  successive  elements  in  the  array.  In  most  cases,  the  bandwidth  of  s(t)  and  the  di¬ 
mensions  of  the  array  are  small  enough  so  that  this  condition  is  satisfied.  Therefore,  s(t) 
remains  constant  during  the  tmn  (#0,^0)  seconds  and  the  approximation  s(t  +  tmn 
*  s(t)  holds.  In  addition,  the  carrier  term  e,0,c'  can  be  removed  since  the  phase  delay  ex¬ 
perienced  by  the  carrier,  is  of  interest.  The  signals  received  at  the  (m,/7),h 

element  then  becomes 

XJt)  =  sU)e'^MI+nm,t(,)  (32) 


33 


and  can  be  rewritten  in  vector  form  given  by: 


X(t)  =  s(t)v  +  n(t)  (33) 

where  v  is  the  steering  vector  associated  with  the  direction  (#0,^0)  of  the  desired  signal 
and  n(t)  is  noise  vector.  The  steering  vector  y  contains  the  AOA  information  and  is 
represented  as  a  one-dimensional  MNx  1  vector  given  by 

y  _  gJmc‘ MliQi’fo)  gi^'ln  giMc'Mn^A) >  gje>c'\N (^>-4) )  ^3 4^ 

where  r  represents  the  matrix  transposition.  The  noise  vector  n(t)  is  written  as  a  one  di¬ 
mensional  MNx  1  vector  given  by 

n(t)  =  [nu(t\  nmX(t ),  ... ,  nm(t)9  njt\  ... ,  nxln(t)9  nw(t)9  ...  ,  nUN(t)f.  (35) 

The  beamformer  output  T(/)can  be  written  as 

Y(t)  =  wX(t)  (36) 

where  *  represents  the  complex  conjugate  operation,  w  is  a  one-dimensional  1  xMN  vec¬ 
tor  given  by 

W*(0  =  [»,(*)>  wmX(t)9  ... ,  wm(t)9  w\H{t)9  ... ,  wMn(t)9  wXN(t)9  ...  ,  (37) 


and  A'(Z)  is  an  MNx  1  vector  given  by 


m  =  [Xn{t\  Xmx(t\  ... ,  Xm{t\  Xxm{t)9  ... ,  xw(t)9  ... ,  XMiV(0f.  (38) 


The  complex  weights  w  are  chosen  in  order  to  create  a  maximum  array  factor  at 
0  —  6o  and  f=<fo  given  by  Equation  (23).  Hence,  the  narrow-band  beamformer  in  Fig. 

1 8  provides  a  means  to  create  a  maximum  array  factor  in  the  desired  direction  by  select¬ 
ing  the  appropriate  complex  weights.  This  works  well  when  the  AOA  is  fixed.  In  cases 
when  the  point  source  is  mobile  (e.g.,  an  UAV),  the  AOA  changes  and  an  adaptive  beam- 
former  is  required,  which  is  presented  in  the  next  section. 


E.  ADAPTIVE  BEAMFORMING 

The  AOA  changes  when  the  point  source  is  moving  and  an  adaptive  beamformer 
is  required  to  steer  the  main  lobe  towards  the  point  source  when  the  latter  is  mobile. 
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Various  adaptive  algorithms  exist,  such  as  the  least  mean  squares  (LMS),  recur¬ 
sive  least  squares  (RLS)  and  constant  modulus  algorithm  (CMA).  Both  the  LMS  and 
RLS  algorithms  require  knowledge  of  the  desired  signal  s(t)  (or  something  closely  corre¬ 
lated  to  it),  while  the  CMA  is  a  blind  adaptive  algorithm,  which  requires  no  knowledge  of 
the  desired  signal.  The  LMS  algorithm  is  based  on  the  steepest-descent  optimization 
method  that  recursively  computes  and  updates  the  complex  weights.  It  is  the  most  com¬ 
mon  adaptive  algorithm  for  continuous  adaptation  and  has  been  well  studied  and  is  well 
understood  [8].  In  this  thesis,  the  LMS  algorithm  is  selected  for  the  adaptive  beamform¬ 
ing  process  due  to  its  simplicity. 

Figure  19  shows  the  block  diagram  of  a  narrow-band  adaptive  beamformer  using 
the  LMS  algorithm.  It  can  be  used  to  steer  the  main  beam  towards  the  point  source  while 
simultaneously  suppressing  interfering  signals  through  adjustment  of  the  complex 
weights. 


Least  Mean  Squares  Beamformer 

(LMS)  Algorithm 


Figure  19.  MxN  narrow-band  adaptive  beamformer. 

The  error  e(t)  is  defined  to  be: 

e(t)  =  r{t)-Y(t)  (39) 
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where  r(t)  is  the  reference  signal  and  Y(t)  is  the  output  of  the  beamformer.  In  practice, 
r(t)  needs  to  be  closely  correlated  to  the  desired  signal  s(t). 

The  complex  weights  are  adjusted  and  applied  to  the  signal  received  by  the  array 
elements  in  an  iterative  manner.  For  each  iteration,  the  complex  weights  are  chosen  to 
minimize  the  mean-square  error  between  the  beamformer  output  T(/)  and  the  reference 
signal  r(t)  given  by  [8] 

£[^(0]  =  £[(r(/)-y(0)2]  (40) 

which  can  be  expanded  to  the  form  given  by 

E  \f  (/)]  =  E  [r2  (0]  -  2w"  E[r(t)X(t)]  +  w"  E[X(t)X"  (/)]  w,  (41) 

where  £[•]  is  the  expectation  operator.  The  minimum  mean-square  error  (MMSE)  is 
obtained  by  setting  the  gradient  vector  of  Equation  (41)  with  respect  to  w  to  zero,  giving 

y„£ [V(0]  =  -2 E[r(t)X(t)]  +  2E[xm"{t)\  w.  (42) 

According  to  the  method  of  steepest  descent,  the  value  of  the  weight  vector  at  time  (/  + 1) 
is  updated  as  follows 

w(l  +  \)  =  ,v(0  +  ^//[-Y,£[^(0]]  (43) 

where  JJ  is  the  convergence  factor.  By  substituting  from  Equation  (42)  into 

the  above  equation,  we  have 

H it  + 1)  =  M<0  + 1//  [2  £  [r(0X(0]  -  2  E[X(t)X"  (0]  *] 

=  w{t)  +  n[r{t)X{t)-X{t)XH{t)w\ 

which  can  be  simplified  to  become: 

w(t  + 1)  =  w(t)  +  jjX(t)e\t)  (44) 
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where  *  represents  complex  conjugation  and  the  instantaneous  estimates  given  by 
[r(/)A'(/)]  and  [A’(/)Ar//(/)]  have  been  used  instead  of  E[r(t)X(t)]  and 

E[X(t)XH(t)^,  respectively.  The  term  E^X(t)XH (t)~^  is  commonly  known  as  the  ar¬ 
ray  correlation  matrix  R  and  denotes  the  correlation  between  the  signals  received  by  the 
various  array  elements.  The  convergence  factor  jj  is  related  to  the  eigenvalues  of  XX" . 

If  the  eigenvalues  of  XX"  are  widely  spread,  convergence  will  be  slow.  For  convergence 
of  the  LMS  algorithm,  the  value  of  ju  needs  to  be  in  the  range  given  by  [20] 

0</7  <2/  Amax  (45) 

where  Aniax  is  the  maximum  eigenvalue  ofXX". 

The  narrow-band  adaptive  beamformer  presented  can  be  used  to  beamform  and 
steer  the  main  beam  of  the  array  towards  the  mobile  point  source.  However,  knowledge 
about  the  AOA  of  the  desired  signal  is  a  prerequisite  for  the  beamforming  process.  As 
such,  methods  of  determining  the  AOA  of  the  desired  signal  are  presented  in  the  next  sec¬ 
tion. 

F.  AOA  ESTIMATION  AND  TRACKING 

The  knowledge  about  the  AOA  of  the  desired  signal  is  a  prerequisite  for  the 
beamforming  process.  Although  it  is  not  the  focus  of  this  thesis  to  provide  an  in-depth 
study  of  existing  AOA  estimation  and  tracking  methods,  they  are  briefly  introduced  in 
this  section  for  completeness.  The  intention  is  to  provide  an  overview  of  the  common 
AOA  estimation  techniques  as  summarized  by  [17].  Methods  to  track  the  estimated  AOA 
are  also  discussed. 

1.  Spectral  Estimation  Methods 

Spectral  estimation  methods  estimate  the  AOA  by  analyzing  the  spatial  spectrum 
and  then  determining  the  local  maxima.  According  to  [17],  one  of  the  earliest  methods  is 
the  Barlett  method  which  applies  equal  weighting  on  each  array  element  and  estimates 
the  mean  input  power.  The  direction  that  gives  the  highest  mean  power  is  the  AOA  of 
the  desired  signal  source.  The  antenna  beamwidth  needs  to  be  very  narrow  and  the 
sidelobe  levels  need  to  be  low  with  respect  to  the  main  lobe  in  order  for  this  method  to 
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work  well.  This  is  because  the  input  power  is  contributed  from  both  the  main  lobe  and 
the  sidelobes  of  the  beam.  A  high  sidelobe  level  may  mislead  the  algorithm  in  deducing 
that  the  AOA  of  the  source  is  in  the  direction  in  which  the  sidelobe  is  pointing  and  not  the 
direction  in  which  the  main  lobe  is  pointing. 

Another  spectral  analysis  method  is  the  minimum-variance  distortionless  response 
(MVDR),  which  finds  the  maximum  likelihood  estimate  of  the  power  arriving  from  a 
point  source  in  direction  0O  assuming  all  other  sources  as  interferences.  The  array 
weights  are  adjusted  to  cancel  the  interference  signals  and  maintain  a  unity  power  gain  at 

2.  Linear  Prediction  Method 

This  method  estimates  the  output  of  one  sensor  using  linear  combinations  of  the 
remaining  sensor  outputs  and  minimizes  the  mean  square  prediction  error,  which  is  the 
error  between  the  estimated  and  actual  output.  It  obtains  the  array  weights  by  minimizing 
the  mean  output  power  of  the  array  subject  to  the  constraint  that  the  weight  on  the  se¬ 
lected  array  element  is  unity. 

3.  Multiple  Signal  Classification  (MUSIC) 

MUSIC  is  a  type  of  eigenstructure  method  that  relies  on  the  properties  of  the  array 
correlation  matrix  R.  The  space  spanned  by  the  eigenvectors  of  R  may  be  partitioned  into 
a  signal  subspace  and  a  noise  subspace  [17].  The  steering  vectors  corresponding  to  the 
directional  sources  are  orthogonal  to  the  noise  subspace  and  are  found  within  the  signal 
subspace.  Several  variations  of  MUSIC  exists,  but  in  its  standard  form,  the  method  esti¬ 
mates  the  noise  subspace  from  the  samples  collected  at  the  array  elements.  Once  the 
noise  subspace  has  been  estimated,  a  search  for  AOA  of  the  desired  signal  is  made  by 
looking  for  steering  vectors  that  are  as  orthogonal  to  the  noise  subspace  as  possible. 

4.  Estimation  of  Signal  Parameters  via  Rotational  Invariant  Techniques 
(ESPRIT) 

This  method  requires  two  arrays  with  identical  element  displacement  vectors.  In 
other  words,  the  elements  in  the  first  array  are  displaced  by  the  same  distance  in  the  same 
direction  from  each  other  as  the  elements  in  the  second  array.  These  two  arrays  can  be 
two  separate  arrays  or  part  of  a  larger  array.  Measurements  are  taken  from  the  two  arrays 
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and  their  array  correlation  matrices  are  estimated.  The  eigenvalues  and  eigenvectors  are 
then  computed,  and  the  directionality  of  the  source  is  determined. 

5.  Methods  of  T  racking  AO  A 

When  the  signal  source  is  mobile  (e.g.,  an  UAV),  the  AOA  of  the  desired  signal 
source  changes,  and  the  next  task  is  to  track  the  change  in  the  AOA  and  steer  the  beam 
towards  the  mobile  signal  source.  Several  approaches  have  been  proposed  to  track  the 
AOA  of  moving  targets  using  passive  arrays  [21-24].  Passive  arrays  do  not  initiate  any 
signal  transmission  and  are  always  in  the  “receiving”  mode.  A  general  assumption  used 
in  these  techniques  is  that  the  targets  can  be  approximated  as  stationary  during  a  limited 
integration  time  interval,  and  the  MUSIC  algorithm  is  applied  to  estimate  the  AOA  pa¬ 
rameters  at  each  time  interval.  Another  approach  proposed  by  [21]  uses  a  maximum  like¬ 
lihood  algorithm  while  that  proposed  by  [22]  uses  a  recursive  tracking  algorithm  that  uses 
the  AOA  estimates  from  the  most  recent  array  data  to  update  the  existing  one.  The  ap¬ 
proach  proposed  by  [23]  attempts  to  minimize  the  distance  between  the  components  of 
the  true  and  the  estimated  array  correlation  matrices  to  provide  an  estimate  of  the  AOA. 

In  summary,  some  common  AOA  estimation  and  tracking  methods  were  dis¬ 
cussed  briefly.  However,  the  methods  introduced  in  this  section  address  the  signal  proc¬ 
essing  aspects  of  determining  the  AOA  only  and  there  is  also  a  need  to  address  the  ques¬ 
tion  of  what  type  of  antenna  beam  is  appropriate  for  AOA  estimation  and  tracking  pur¬ 
poses  and  whether  active  or  passive  detection  is  suitable.  These  questions  are  addressed 
in  the  next  section. 

G.  DETECTION  AND  TRACKING  OF  UAV 

After  presenting  the  methods  to  estimate  and  track  the  AOA  of  the  desired  signal 
sent  by  the  UAV,  this  section  discusses  some  of  the  considerations  when  using  arrays  for 
detection  and  tracking  of  the  UAV. 

1.  Active  and  Passive  Detection 

Figure  20  shows  the  difference  between  active  and  passive  detection.  For  active 
detection,  a  sensor  cluster  constantly  searches  for  the  presence  of  the  UAV  by  sending  a 
predetermined  signal  and  measuring  the  returned  signal  characteristics.  When  the  UAV 
has  been  detected  and  its  position  determined,  the  sensor  cluster  initiates  beamforming 
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and  tracking  of  the  UAV  and  performs  the  task  of  transmitting  the  collected  data  to  the 
UAV  through  wireless  means. 

Sensor  nodes  typically  have  very  limited  energy  available  since  they  are  battery- 
operated.  From  an  energy  dissipation  perspective,  the  passive  detection  approach  is  pre¬ 
ferred.  A  passive  sensor  cluster  constantly  searches  for  a  predefined  signal  sent  by  the 
UAV.  This  predefined  reference  signal  may  be  sent  continuously  (or  at  regular  time  in¬ 
tervals)  towards  the  sensor  cluster.  The  sensor  cluster  then  uses  the  received  signal  to  de¬ 
termine  the  AOA  of  the  signal  source  and  initiates  beamforming  and  tracking  of  the 
UAV. 

The  advantage  of  employing  passive  detection  is  that  the  sensor  cluster  is  not  re¬ 
quired  to  transmit  signals  for  detection  of  the  UAV,  hence  passive  detection  by  the  sensor 
cluster  is  adopted  in  this  thesis. 


UAV 

aar-_'-===i6 


03T-C- 


Active  Detection 


Passive  Detection 


Figure  20.  Active  and  passive  detection  of  UAV  using  a  sensor  cluster. 


2.  Detection  and  Tracking  of  UAV 

A  broad  beam  is  especially  useful  for  detection  of  UAV  because  of  the  larger  de¬ 
tection  region  being  created  as  illustrated  in  Fig.  21 .  The  detection  region  becomes  larger 
as  the  altitude  increases,  making  the  broad  antenna  beam  a  very  suitable  candidate  for 
UAV  detection.  Therefore,  it  is  appropriate  to  use  as  few  array  elements  as  possible  to 
generate  a  broad  beam  to  detect  the  UAV  in  this  case. 
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Detection  region  lor  narrow  beam  u  th  f>. 

Detection  region  lor  broad  beam  with  (>%,*  (• 

Figure  21 .  Generating  a  broad  antenna  beam  to  increase  the  size  of  the  UAV  detection 

region. 

It  is  often  not  desirable  to  use  a  broad  beam  to  track  the  UAV  because  of  the 
lower  tracking  resolution  it  provides  compared  to  that  of  a  narrower  beam.  Figure  22  il¬ 
lustrates  how  the  tracking  resolution  of  a  broader  beam  degrades  as  the  altitude  decreases. 
As  shown  in  the  figure,  it  is  not  possible  to  resolve  or  differentiate  between  the  two 
UAVs  at  altitude  /z?  because  both  of  their  positions  reside  within  the  3-dB  beamwdith  of 
the  broader  beam.  Hence,  it  is  more  appropriate  in  this  case  to  use  as  many  array  ele¬ 
ments  as  possible  in  order  to  generate  a  narrow  beam  to  provide  the  tracking  resolution 
required. 
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UAV 


Figure  22.  Generating  a  narrow  beam  to  increase  the  tracking  resolution  of  the  UAV. 

H.  SUMMARY 

This  chapter  presented  the  fundamental  concepts  required  for  each  sensor  cluster 
to  form  an  array-like  transmission  to  the  UAV.  In  particular,  the  array  factor  and  power 
gain  equations  of  antenna  array  were  presented,  followed  by  a  literature  review  on  some 
findings  on  random  arrays.  This  was  followed  by  the  discussion  on  the  narrowband 
adaptive  beamformer  and  the  various  AOA  estimation  and  tracking  algorithms  using 
passive  arrays.  The  significance  of  the  antenna  beamwidth  for  UAV  detection  and  track¬ 
ing  was  also  highlighted. 

Up  to  this  point,  the  relevant  discussion  has  been  presented  in  order  to  develop  a 
simulation  model  based  on  the  operation  scenario  shown  in  Figures  1  and  2.  The  simula¬ 
tion  model  and  simulation  results  are  presented  in  the  next  chapter.  The  objective  is  to 
study  the  use  of  sensor  clusters  for  adaptive  beamforming  through  simulation  and  com¬ 
pare  the  simulation  results  with  the  theoretical  results  discussed  in  Section  C  of  Chapter 
III. 
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IV.  SIMULATION  RESULTS 


A  sensor  cluster  differs  from  a  conventional  array  in  the  following  aspects:  1)  an 
unknown  number  of  sensor  nodes  used  for  beamforming,  2)  random  positioning  of  the 
sensor  nodes,  3)  position  errors  due  to  mobility  of  sensor  nodes  and  4)  frequent  failures 
of  sensor  nodes.  These  are  issues  found  in  the  design  and  practical  implementation  of  a 
wireless  sensor  network. 

This  chapter  aims  to  study  the  effects  of  the  above  mentioned  issues  on  the  an¬ 
tenna  beam  used  for  data  transmission  from  a  sensor  cluster  to  an  UAV.  First,  the  simu¬ 
lation  scenario  given  in  Figures  1  and  2  is  presented  in  greater  detail.  This  is  followed  by 
a  discussion  on  the  simulation  model  developed  here  to  study  the  use  of  sensor  cluster  for 
adaptive  beamforming.  Next,  simulation  results  are  presented  and  compared  to  the  theo¬ 
retical  findings  of  random  arrays  discussed  in  Section  C  of  Chapter  III.  The  effect  of 
UAV  speed  on  the  convergence  of  the  adaptive  beamforming  algorithm  is  also  discussed. 

A.  SIMULATION  SCENARIO 

This  section  describes  the  simulation  scenario  presented  in  Figures  1  and  2.  The 
three-dimensional  and  its  equivalent  two-dimensional  simulation  scenario  are  shown  in 
Figures  23  and  24,  respectively.  Both  figures  show  the  flight  path  of  an  UAV  which  at¬ 
tempts  to  receive  the  data  collected  by  the  sensor  nodes  as  it  flies  over  them. 

The  sequence  of  events  for  this  simulation  scenario  is  described  as  follows.  After 
the  sensor  nodes  are  deployed  randomly  over  an  area  of  interest,  they  are  grouped  to¬ 
gether  into  sensor  clusters  and  a  primary  node  is  elected  for  each  cluster  based  on  meth¬ 
ods  described  in  Section  B  of  Chapter  II.  It  is  further  assumed  that  the  positions  of  the 
secondary  nodes  within  the  sensor  cluster  are  determined  using  methods  presented  in 
Section  C  of  Chapter  II.  Next,  an  UAV  flies  over  the  sensor  field  and  sends  a  predeter¬ 
mined  desired  signal  to  these  sensor  clusters  continuously.  Assuming  that  each  sensor 
cluster  has  knowledge  of  the  AOA  of  the  desired  and  interfering  signals  (based  on  meth¬ 
ods  discussed  in  Chapter  III,  Section  F),  it  adapts  to  the  change  in  the  AOA  of  desired 
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z-axis  UAV 


Figure  23.  Three-dimensional  simulation  scenario.  The  UAV  sends  a  predetermined  desired 
signal  towards  the  sensor  clusters  in  the  direction  of  00  =  +20  to  -20  at  ft  =  +45 .  The 

two  interference  signals  arrive  at  the  sensor  clusters  at(#  =  ±60  ,(f  =  45  ). 


Array  Normal 


Figure  24.  Two-dimensional  diagram  showing  the  scenario  in  Fig.  23  taken  at  the  plane  AA\ 
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signal  sent  from  the  UAV  and  steers  the  antenna  beam  towards  the  UAV  using  the  LMS 
algorithm.  The  following  assumptions  apply  to  this  simulation  scenario: 


1 .  Passive  detection  of  the  UAV  by  the  sensor  cluster  array  as  discussed  in 
Chapter  III,  Section  G.  1 . 

2.  The  desired  complex  modulating  signal  s(t)  sent  by  the  UAV  is  known  to 
the  sensor  cluster  array  and  is  the  same  as  the  reference  signal  r(t)  used  for 
the  LMS  algorithm. 

3.  Each  primary  node  has  perfect  knowledge  of  its  own  position  and  the  posi¬ 
tion  of  the  secondary  nodes  within  the  sensor  cluster  prior  to  beamforming. 

4.  The  AO  A  of  the  desired  signal  (#0,^0 )  is  known  to  the  sensor  cluster.  The 
elevation  angle  #0  spans  from  +20  to  -  20  and  the  azimuth  angle  j){)  is 
fixed  at  +45\ 

5.  The  AO  A  of  the  first  interference  signal  (#,  ,^, )  is  known  to  the  sensor  clus¬ 
ter  and  is  fixed  at  elevation  angle#,  =  +60  and  azimuth  angle <f{  =  +45  . 

6.  The  AOA  of  the  second  interference  signal  (#,,^2)  is  known  to  the  sensor 
cluster  and  is  fixed  at  elevation  angle#,  =  -60  and  azimuth  angle </>2  =  +45  . 

7.  There  is  perfect  frequency,  phase  and  data  synchronization  among  all  the 
sensor  nodes  within  the  sensor  cluster  and  between  the  UAV  and  each  sen¬ 
sor  cluster. 

8.  There  is  perfect  coordination  among  the  primary  node  and  the  secondary 
nodes  within  a  sensor  cluster  for  beamforming. 

9.  All  sensor  cluster  arrays  are  designated  to  operate  at  a  carrier  frequency  of 
500  MHz  (i.e.,  X  =  0.6  m). 

10.  All  sensor  nodes  are  modeled  as  isotropic  antenna  elements  and  mutual¬ 
coupling  between  sensor  nodes  is  not  considered. 
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The  description  of  the  simulation  scenario  and  its  assumptions  set  the  baseline  for 
the  development  of  the  simulation  model  presented  in  the  next  section. 

B.  SIMULATION  MODEL 

Figure  25  shows  the  simulation  model  used  in  this  thesis.  The  model  was  imple¬ 
mented  in  MATLAB  programming  language  and  the  program  descriptions  are  found  in 
the  Appendix. 


XuO  xnJt)  xm(t)  xln(t)  xMn(t)  xIN(t)  -W V 


Least  Mean  Squares  Beam  former 

(LMS)  Algorithm 


Figure  25.  Simulation  model  for  MxN  adaptive  narrow-band  beamformer.  The  complex 
weights  are  adjusted  by  the  LMS  algorithm  and  its  outputs  are  applied  to  the  array 

elements  to  form  a  beam. 

As  mentioned  in  the  previous  section,  the  desired  complex  modulating  signal  s(t) 
sent  by  the  UAV  is  known  to  the  sensor  cluster  array  and  becomes  the  reference  signal 
for  the  LMS  algorithm.  The  array  signal  generator  generates  the  array  outputs  X\t) 
given  by 

X\t)  =  s(t)V  +  n\t)  +  z,  {t)u[  +  z2  {t)u[  (46) 

which  is  a  modified  version  of  Equation  (33).  For  convenience,  the  desired  complex 
modulating  signal  s(t)  is  assumed  to  be  a  real-valued  modulating  signal  given  by 
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,  .  2k ft 

.9(0  =  Sin  ——+y/a 

J  s 


(47) 


where/ is  the  frequency  in  Hz,  t  is  the  discrete  time  in  seconds,  f  is  the  number  of  sam¬ 
ples  per  second  (or  sampling  rate)  and  y/o  is  a  uniform  random  variable  in  the  range  0  to 
271.  The  steering  vector  y'  associated  with  the  direction  of  the  desired  signal  is  a 

two  dimensional  MxN  vector  given  by 


...  ...  gJMA/hWaJo) 


gJ°*  W'WiA))  ...  gj**MNWQJ0)  ... 


(48) 


which  is  a  modified  version  of  Equation  (34)  in  which  v  was  an  MNx  1  vector.  The  one¬ 
dimensional  representation  is  a  more  convenient  form  of  mathematical  representation 
while  the  two-dimensional  representation  is  preferred  for  MATLAB  implementation  and 
is  adopted  in  this  chapter.  Similarly,  n(t)  is  an  MxN  noise  matrix  given  by 

«n(0  •••  «.,(0  ••• 

n\t)=  njt)  ...  «.(/)  -  «,,„(/)  (49) 

_«,,v(0  •••  ••• 

which  is  a  modified  version  of  Equation  (35). 

The  two  interfering  modulating  signals z,(/)  and z2(t)  with  the  corresponding  steer¬ 
ing  vectors  u[  and  ur2,  respectively,  are  also  included  in  Equation  (46).  The  objective  is 
to  steer  the  main  beam  in  the  direction  of  the  desired  signal  and  suppress  the  interference 
signals  by  creating  nulls  in  their  directions.  The  first  interfering  modulating  signal  z,(/) 
is  given  by 
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,,  .  2a  ft 

z,(0  =  sin|  — —  +  <-/, 

Js 


(50) 


where  f/j  is  an  uniform  random  variable  in  the  range  0  to  2rc.  The  steering  vector  u[  as¬ 
sociated  with  the  direction  (dr/ )  of  the  first  interference  signal  is  an  MxN  matrix  given 


J<# 1 1  (^  ) 


Ja*M  1  -d  > 


,MibV  > 


The  second  interfering  modulating  signal  z2(t)  is  given  by 


,  .  .  2k  ft 

z,(0  =  sin  — —  +  y/2 

J  s 


where  j/2  is  a  uniform  random  variable  in  the  range  0  to  2n.  The  steering  vector  u[  asso¬ 
ciated  with  the  direction  )  of  the  second  interference  signal  is  an  MxN  matrix  given 


»/**!  \(°2’*2) 


J°*M  |  <^2  -^2  > 


,jotliN(&2^) 


J**mN  <?2  -&  > 


The  array  output  X\t)  is  an  MxN  matrix  given  by 


*.,(')  -  *.,(0  -  **,(') 


X'(t)=  Xjt)  -  xjt)  ...  x,,„(0 


*„(0  -  ^„,(o  -.  xm(t) 


which  is  a  modified  version  of  Equation  (38). 
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The  output  vector  of  the  beamformer  Y'(t)  is  given  by 


Y\t)  =  W\t)*X\t) 

where  •  denotes  dot  product  and  w'(t)  is  an  MxN weight  matrix  given  by 


(55) 


w'‘(0  = 


%(')  •••  4,(0  •••  4,(0 


4(0  •••  4,(0  4,(0 

4(0  •••  4(0  •••  4,(0 


(56) 


which  is  a  modified  version  of  Equation  (37).  Therefore,  Y\t)  becomes  an  an  MxN  ma¬ 
trix  given  by 

4,(0  -  T,(0  -  Ym(t) 


no = 


yj o  -  uo  -  4( o 


4(0  -  4(0  -  4(0 


(57) 


The  beamformer  output  Y(t)  given  in  Equation  (36)  is  obtained  by  summing  all  the  ma¬ 
trix  elements  within  the  Y\t)  matrix  given  in  Equation  (57). 

The  simulation  procedure  is  to  sample  s(t),  z/(t)  and  Z2(t)  at  fy  and  the  complex 
weights  are  computed  for  each  iteration  t  using  the  LMS  algorithm.  The  algorithm  is  said 
to  have  converged  when  £^e2(/)J  <  0.1,  which  is  defined  as  the  convergence  threshold. 

The  complex  weights  generated  are  then  used  to  compute  the  array  factor  using  Equation 
(23)  and  plot  the  beam  pattern  (or  power  gain)  using  Equation  (24). 


C.  EFFECTS  OF  SENSOR  NODE  DENSITY 

In  this  section,  the  effects  of  the  number  of  sensor  nodes  on  the  antenna  beam 
generated  by  a  sensor  cluster  are  presented.  The  location  of  each  node  is  assumed  to  be 
uniformly  distributed  within  a  1 .2x1 .2  nr  plane;  hence,  changing  the  number  of  sensor 
nodes  within  the  sensor  cluster  is  equivalent  to  changing  the  density  of  sensor  nodes 
within  the  cluster.  The  simulation  was  performed  by  assuming  that  the  AOA  of  the  de- 
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sired  signal  is  at  6{)  =  0  and  f0  =  +45  .  One  hundred  sets  of  sensor  clusters  are  generated 
for  each  sensor  node  density  and  the  corresponding  power  gains  computed.  The  average 
power  gain  is  then  obtained  by  averaging  these  results.  This  method  of  obtaining  the  av¬ 
erage  power  gain  is  adopted  throughout  the  chapter. 

The  interference  rejection  level  Ir  indicates  how  well  the  antenna  beam  can  sup¬ 
press  the  interference  level.  A  high  Ir  level  is  desirable  in  the  direction  of  the  interfer¬ 
ence  signal.  Figure  26  shows  the  plots  of  the  antenna  beam  pattern  for  sensor  cluster 
densities  ranging  from  3x3  to  16x16  sensor  nodes.  It  can  be  seen  that  the  3-dB  beam- 
width  and  the  shape  of  the  main  lobe  remain  relatively  unchanged  as  the  number  of  sen¬ 
sor  nodes  change.  However,  the  maximum  average  power  gain  Gmax  and  the  interference 
rejection  level  Ir  decreases  as  the  number  of  sensor  nodes  is  decreased  to  less  than  7x7. 
Figure  27  shows  that  Gmax  has  reduced  from  0  dB  for  the  7x7  sensor  cluster  to  -3.3  dB 
for  the  3x3  sensor  cluster  while  Fig.  28  shows  that  the  interference  rejection  is  reduced 
from  55  dB  for  the  9x9  sensor  cluster  to  1 7  dB  for  the  3x3  sensor  cluster. 

These  results  are  consistent  with  those  discussed  in  Chapter  III,  Section  C.l  where 
the  results  derived  by  [7]  suggested  that  the  average  sidelobe  level  can  be  reduced  by  in¬ 
creasing  the  density  of  the  sensor  nodes  within  the  cluster.  In  our  case  here,  increasing 
the  number  of  sensor  nodes  enabled  the  reduction  of  the  average  sidelobe  levels  and 
hence  suppressed  the  interference  signals  more  effectively. 
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Figure  26.  Average  power  gain  of  the  beams  generated  by  16x16,  7x7,  6x6,  5x5,  4x4  and 
3x3  sensor  clusters  at  (0O,^O)  =  (0  ,45  ) :  (a)  Polar  plots  and  (b)  X-Y  plots  showing  the 

decrease  in  maximum  average  power  gain  and  interference  rejection  level  as  the  number 

of  sensor  nodes  decreased. 
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Figure  27.  The  maximum  average  power  gain  as  a  function  of  the  number  of  sensor  nodes. 


Figure  28.  The  interference  rejection  level  as  a  function  of  the  number  of  sensor  nodes. 
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D.  EFFECTS  OF  POSITION  ERRORS 

It  was  assumed  that  the  primary  node  has  perfect  information  of  the  initial  loca¬ 
tions  of  the  secondary  nodes  (using  the  location  discovery  methods  described  in  Chapter 
II)  prior  to  beamforming.  In  reality,  there  may  be  errors  associated  with  these  location 
discovery  methods.  Moreover,  the  positions  of  these  mobile  sensor  nodes  may  change 
without  the  knowledge  of  the  primary  node  during  the  beamforming  process.  These 
cause  position  errors  (or  phase  errors).  The  consequence  is  that  the  complex  weights 
used  for  beamforming  are  being  adjusted  without  accounting  for  the  change  in  the  loca¬ 
tions  of  the  sensor  nodes.  This  affects  the  average  antenna  array  factor  and,  hence,  the 
average  beam  pattern  (or  average  power  gain). 

In  this  section,  the  effects  of  position  errors  on  the  antenna  beam  generated  by  a 
7x7  sensor  cluster  are  presented.  It  is  assumed  that  the  sensor  nodes  move  much  slower 
than  the  speed  of  the  UAV  and  move  within  a  distance  of  ±A d  from  their  initial  loca¬ 
tions  known  to  the  primary  node.  Therefore,  Ad  is  regarded  as  the  position  error.  In 
particular,  five  values  of  Ad  are  considered:  0.1A,  0.2 A,  0.3 A,  0.4A,  and  0.5A,  where 
/i  =  0.6  m.  The  simulation  was  performed  by  assuming  that  the  AOA  of  the  desired  sig¬ 
nal  is  at0o  =  0  and  (j>()  =  +45°. 

Figure  29  shows  the  plots  of  the  antenna  beam  pattern  for  different  values  of  Ad. 
It  shows  that  the  maximum  average  power  gain  and  the  3-dB  beamwidth  are  not  affected 
significantly  by  the  magnitude  of  Ad.  This  result  is  better  than  that  discussed  in  Chapter 
III,  Section  2,  where  it  was  shown  that  a  phase  error  of  ±40  is  required  to  cause  a  3-dB 
reduction  in  the  maximum  average  power  gain.  In  this  case,  the  maximum  average 
power  gain  remained  the  same  even  though  the  phase  error  was  0.5A  (or  ±90  ). 

It  can  be  seen  from  Fig.  29  that  the  interference  rejection  level  Ir  decreased  when 
Ad  increased.  Figure  30  shows  that  the  value  of  Ir  decreased  from  28  dB  when  A d  = 

0.1  A  to  23  dB  when  A d  =  0.2A,  a  reduction  of  5  dB  when  Ad  is  doubled.  Similarly,  the 
value  of  Ir  reduced  by  5  dB  from  23  dB  to  17  dB  when  Ad  increased  from0.2A  to  0.4/1. 
This  is  1  dB  less  than  that  discussed  in  Chapter  III,  Section  2. 
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- No  position  errors  - Ad-  OJA  (1 8  cm)  2-D  plots  generated  at  4>  =  4.< 

- Ad =  0. 1A  (6  cm)  Ad-  0.4A  (24  cm) 

Ad-  0.2A  (12  cm)  - Ad-  OJA (30 cm) 

Figure  29.  Average  power  gain  of  the  beams  generated  by  the  7x7  sensor  cluster  in  the  pres¬ 
ence  of  position  errors  A d=  0.1A,  0.2A,  0.3A,  0.4A,  and  0.5A  at  (#<,,$,)  =  (0',45'): 

(a)  Polar  plots  and  (b)  X-Y  plots  showing  the  decrease  in  interference  rejection  level  as 

the  position  errors  increased. 
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E.  EFFECTS  OF  SENSOR  NODE  FAILURES 

Sensor  nodes  fail  frequently  in  wireless  sensor  networks  due  to  depleted  battery 
energy  or  damaged  by  environmental  conditions.  The  effects  of  sensor  node  failures  on 
the  antenna  beam  being  generated  by  the  sensor  cluster  are  presented  in  this  section  using 
a  7x7  sensor  cluster. 

The  simulation  was  performed  by  assuming  that  the  AOA  of  the  desired  signal  is 
at  00  =  0  and  (/>0  =  +45°.  Figure  3 1  shows  the  antenna  beam  pattern  produced  by  the  7x7 

sensor  cluster  for  five  cases:  1)  no  element  fails,  2)  five  elements  fail,  3)  ten  elements 
fail,  4)  fifteen  elements  fail  and  5)  twenty  elements  fail.  It  shows  that  the  3-dB  beam- 
width  and  the  shape  of  the  main  lobe  remain  relatively  unchanged  in  the  presence  of  sen¬ 
sor  node  failures.  However,  there  is  some  reduction  in  the  maximum  average  power  gain 
as  the  number  of  sensor  node  failures  increases.  For  example,  the  maximum  average 
power  gain  is  reduced  by  4.65  dB  when  twenty  sensor  nodes  failed  compared  to  the  case 
when  no  sensor  nodes  failed.  This  is  expected  since  the  magnitude  (of  the  complex 
weights)  contribution  from  the  sensor  nodes  decreases  due  to  sensor  node  failures.  Addi¬ 
tionally,  the  interference  rejection  levels  decreased  as  the  number  of  sensor  node  failures 
increased,  thus  reducing  the  sensor  cluster’s  capability  to  create  nulls  at  the  AOA  of  the 
two  interfering  signals  at  6  =  ±60  .  The  interference  rejection  level  was  about  52  dB 
when  no  sensor  nodes  failed  and  reduced  to  about  23.75  dB  when  twenty  sensor  nodes 
failed.  Hence,  sensor  node  failures  do  not  affect  the  3-dB  beamwidth  and  the  shape  of 
the  main  lobe  significantly,  but  the  maximum  average  power  gain  and  the  interference  re¬ 
jection  levels  of  the  antenna  beam  generated  by  the  sensor  cluster  are  reduced  considera¬ 
bly. 

Figure  32  shows  the  reduction  in  maximum  average  power  gain  G™*  as  the  num¬ 
ber  of  sensor  node  failures  increases.  The  simulated  results  showed  good  agreement  with 
that  calculated  using  Equation  (28).  It  has  reduced  from  0  dB  in  the  absence  of  sensor 
node  failures  to  about  4.6  dB  when  twenty  sensor  nodes  failed.  Figure  33  shows  that  the 
interference  rejection  has  reduced  from  50  dB  in  the  absence  of  sensor  node  failures  to 
about  24  dB  when  twenty  sensor  nodes  failed. 
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- No  sensor  nodes  failed  1 5  sensor  nodes  failed  2-D  plots  generated  at  =  45 
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Figure  3 1 .  Average  power  gain  of  the  beams  generated  by  the  7x7  sensor  cluster  as  the  num¬ 

ber  of  sensor  node  failures  increased  from  five  to  twenty  at  (0O,^O)  =  (0\45  ) :  (a)  Polar 

plots  and  (b)  X-Y  plots  showing  the  decrease  in  maximum  average  power  gain  and  inter¬ 
ference  rejection  level  as  the  number  of  sensor  node  failures  increased. 
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Figure  32.  The  maximum  average  power  gain  as  a  function  of  sensor  node  failures. 
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Figure  33.  The  interference  rejection  level  as  a  function  of  sensor  node  failures. 
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F.  ADAPTIVE  BEAMFORMING  USING  SENSOR  CLUSTERS 

Adaptive  beamforming  is  implemented  in  simulation  to  steer  the  antenna  beam 
towards  the  UAV  as  it  flies  from  0{]  =  +20  to -20  at  f0  =  +45  using  the  LMS  algorithm. 
This  section  contains  a  series  of  plots  illustrating  the  effects  of  position  errors  and  sensor 
node  failures  on  the  antenna  beam  pattern  (or  average  power  gain)  generated  by  the  7x7 
sensor  cluster.  In  particular,  five  cases  are  presented: 

1.  AO  A  of  the  desired  signal  at  (0O,^O)  =  (20  ,45  ).  No  position  errors  and 
sensor  node  failures  are  experienced  within  the  sensor  cluster. 

2.  AO  A  of  the  desired  signal  at  (0O,^O)=  (lO  ,45  ).  Position  errors  of  Ad  = 
0.1  A  and  five  sensor  node  failures  are  simulated  within  the  sensor  cluster. 

3.  AO  A  of  the  desired  signal  at  (#0,^>0)  =  (0  ,45  ).  Position  errors  of  Ad  = 
0.2A  and  ten  sensor  node  failures  are  simulated  within  the  sensor  cluster. 

4.  AO  A  of  the  desired  signal  at  (#0,^0)  =  (-10  ,45  ).  Position  errors  of  Ad  = 
0.3A  and  fifteen  sensor  node  failures  are  simulated  within  the  sensor  clus¬ 
ter. 

5.  AO  A  of  the  desired  signal  at  (#0,^0)=  (-20  ,45  ).  Position  errors  of  Ad  = 
0.4A  and  twenty  sensor  node  failures  are  simulated  within  the  sensor  clus¬ 
ter. 

For  each  case,  the  antenna  beam  generated  by  the  sensor  cluster  in  the  presence  of 
position  errors  and  sensor  node  failures  are  shown. 

Figure  34  shows  the  two-dimensional  plots  of  the  antenna  average  power  gain 
versus  the  elevation  angle  6  for  case  1 .  It  can  be  seen  that  the  maximum  average  power 
gain  occurs  at  9  =  20',  and  nulls  are  created  at  6  =  ±60  .  This  shows  that  the  sensor  clus¬ 
ter  can  steer  the  main  lobe  towards  the  AOA  of  the  desired  signal  and  suppress  the  inter¬ 
ference  signals  simultaneously.  The  highest  sidelobe  level  is  about  20  dB  below  the  main 
lobe,  and  the  interference  rejection  level  is  about  52  dB. 
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Figure  35  shows  the  two-dimensional  plots  of  the  antenna  average  power  gain 
versus  the  elevation  angle  0  for  case  2.  It  can  be  seen  that  the  sensor  cluster  managed  to 
steer  the  beam  towards  the  AO  A  of  the  desired  signal.  The  shape  of  the  main  lobe  and 
the  3-dB  beamwidth  also  remain  relatively  unchanged  under  these  conditions.  However, 
the  maximum  average  power  gain  has  reduced  by  1  dB  while  the  interference  rejection 
ratio  has  reduced  to  about  25  dB  at  0  =  -60  and  26  dB  at  6  =  60  . 

Figure  36  shows  the  two-dimensional  plots  of  the  antenna  average  power  gain 
versus  the  elevation  angle  0  for  case  3.  It  can  be  seen  that  the  sensor  cluster  still  man¬ 
aged  to  steer  the  beam  towards  the  AO  A  of  the  desired  signal.  The  shape  of  the  main 
lobe  and  the  3-dB  beamwidth  also  remain  relatively  unchanged  under  these  conditions. 
However,  the  maximum  average  power  gain  has  reduced  by  2  dB  while  the  interference 
rejection  ratio  has  reduced  to  about  22  dB  at  0  =  ±60°. 

Figure  37  shows  the  two-dimensional  plots  of  the  antenna  average  power  gain 
versus  the  elevation  angle  0  for  case  4.  It  can  be  seen  that  the  main  lobe  points  towards 
the  AOA  of  the  desired  signal  and  the  shape  of  the  main  lobe  and  the  3-dB  beamwidth 
also  remain  relatively  unchanged  under  these  conditions.  However,  the  maximum  aver¬ 
age  power  gain  is  reduced  by  3.3  dB  while  the  interference  rejection  ratio  is  reduced  to 
about  21  dB  at  0  =  -60  and  22  dB  at  0  =  60". 

Figure  38  shows  the  two-dimensional  plots  of  the  antenna  average  power  gain 
versus  the  elevation  angle  0  for  case  5.  It  can  be  seen  that  the  main  lobe  points  at 
0  =  -15  instead  of  -20  ,  a  main  lobe  pointing  error  of  5°.  This  shows  that  the  sensor 
cluster  is  unable  to  steer  the  beam  towards  the  AOA  of  the  desired  signal  when  A d  = 
0.4A  wavelength  and  20  sensor  nodes  failed.  The  maximum  average  power  gain  has  re¬ 
duced  by  5.5  dB  while  the  interference  rejection  ratio  has  reduced  to  about  20  dB  at 
0  =  -60  ,  and  20.5  dB  at  0  =  60  . 

Figure  39  shows  the  reduction  in  maximum  average  power  gain  Gmax  as  the  posi¬ 
tion  errors  and  sensor  node  failures  increase.  It  has  reduced  from  0  dB  in  case  1  to  about 
5.5  dB  in  case  5.  Figure  40  shows  that  the  interference  rejection  has  reduced  from  50  dB 
in  case  1  to  13  dB  in  case  5.  Both  figures  show  that  as  the  sensor  cluster  beamforms  to- 
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wards  the  UAV,  the  maximum  average  power  gain  and  the  interference  rejection  levels  of 
the  antenna  beam  decrease  when  the  position  errors  and  sensor  node  failures  increase. 

0 


0  (degrees) 


Average  power  gain,  C?(dB) 

(a) 


0  (degrees) 

(b) 

Figure  34.  Average  power  gain  of  the  beams  generated  by  the  7x7  sensor  cluster  in  the 
absence  of  position  errors  and  sensor  node  failures  at  (#0,$,)  =  (20  ,45  ) :  (a)  Polar  plot 
and  (b)  X-Y  plot  showing  the  maximum  average  power  gain  of  0  dB  and  interference 
rejection  level  of  52  dB  and  55  dB  at  6  =  60  and  6  =  -60  ,  respectively. 
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D 

(degrees) 


- Absence  of  position  errors  and  sensor  node  failures 

-  Presence  of  position  errors  0.1  A.)  and  5  sensor  node  failures 

2-D  plots  generated  at  <J>  =  45' 


Figure  35.  Average  power  gain  of  the  beams  generated  by  the  7x7  sensor  cluster  in  the  pres¬ 
ence  of  position  errors  Ad  =  OAA  and  five  sensor  node  failures  at  (#0,^0)  =  (l0  ,45  ) : 

(a)  Polar  plots  and  (b)  X-Y  plots  showing  the  decrease  in  maximum  average  power  gain 

and  interference  rejection  level. 
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0  (degrees) 


- Absence  of  position  errors  and  sensor  node  failures 

-  Presence  of  position  errors  ( [Ad  =  0.2A)  and  1 0  sensor  node  failures 

2-D  plots  generated  at<j»  =  45" 

Figure  36.  Average  power  gain  of  the  beams  generated  by  the  7x7  sensor  cluster  in  the  pres¬ 
ence  of  position  errors  Ad  =  0.2A  and  ten  sensor  node  failures  at  (#„,$,)  =  (0  ,45  ) :  (a) 

Polar  plots  and  (b)  X-Y  plots  showing  further  decrease  in  maximum  average  power  gain 

and  interference  rejection  level. 
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0  (degrees) 


- Absence  of  position  errors  and  sensor  node  failures 

-  Presence  of  position  errors  (Ad=  G.3A)  and  1 5  sensor  node  failures 

2-D  plots  generated  at  4>  =  45' 

Figure  37.  Average  power  gain  of  the  beams  generated  by  the  7x7  sensor  cluster  in  the  pres¬ 
ence  of  position  errors  Ad  =  0.3A  and  fifteen  sensor  node  failures  at  ( d0 ,<(0)  = 

(-10  ,45  ) :  (a)  Polar  plots  and  (b)  X-Y  plots  showing  the  decrease  in  maximum  average 

power  gain  and  interference  rejection  level. 
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(b) 


- Absence  of  position  errors  and  sensor  node  failures 

-  Presence  of  position  errors  {Ad=  0.4A)  and  20  sensor  node  failures 

2-D  plots  generated  at  <j»  =  45" 

Figure  38.  Average  power  gain  of  the  beams  generated  by  the  7x7  sensor  cluster  in  the  pres¬ 
ence  of  position  errors  A d  =  0.4/1  and  twenty  sensor  node  failures  at  (#„,$,)  = 

(-20  ,45  ) :  (a)  Polar  plots  and  (b)  X-Y  plots  showing  the  decrease  in  maximum  average 
power  gain  and  interference  rejection  level  and  the  presence  of  main  lobe  pointing  error. 
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Case  1  Case  2  Case  3  Case  4  Case  5 


Simulation  cases 

Figure  39.  Reduction  in  maximum  power  gain  with  increase  in  Ad  and  sensor  node  failures. 


Case  1  Case  2  Case  3  Case  4  Case  5 


Simulation  cases 

Figure  40.  Reduction  in  interference  rejection  level  with  increase  in  Ad  and  sensor  node 

failures. 
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G.  EFFECTS  OF  UAV  SPEED  ON  THE  ADAPTIVE  ALGORITHM 

The  speed  of  the  UAV  has  some  influence  over  the  performance  of  the  algorithm 
used  for  adaptive  beamforming.  Figure  41  shows  a  plot  of  the  magnitude  of  the  squared 
error  |e2(/)|  versus  the  number  of  samples  Ns.  This  plot  shows  the  convergence  trend  of 

the  algorithm  as  the  number  of  samples  increases.  The  mean  square  error  is  computed  to 
be  0.0905  after  120  iterations  which  is  below  the  convergence  threshold,  an  indication  of 
convergence.  It  was  further  assumed  that  the  algorithm  required  120  samples/s  in  order 
to  reach  convergence  when  the  speed  of  the  UAV  was  40  m/s. 

If  the  speed  of  the  UAV  were  increased  by  a  factor  of  two  (becomes  80  m/s),  the 
sampling  rate  would  be  reduced  by  the  same  factor  (becomes  60  samples/s).  Figure  42 
shows  the  convergence  of  the  LMS  algorithm  for  this  case.  The  mean  square  error  is 
computed  to  be  0.1697  after  60  iterations  which  is  above  the  convergence  threshold. 
Figure  43  shows  how  the  mean  square  error  deteriorates  when  the  speed  of  the  UAV  in¬ 
creases.  This  problem  may  be  solved  by  increasing  the  sampling  rate  of  the  adaptive  al¬ 
gorithm  when  the  speed  of  the  UAV  increases.  As  shown  in  Fig.  43,  the  LMS  can  still 
converge  at  a  sampling  rate  of  360  samples/s  when  the  UAV  speed  is  120  m/s. 


Figure  41 .  Convergence  of  the  LMS  algorithm  at  UAV  speed  of  40  m/s. 
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Magnitude  of  square  error.  \e?(l)\ 


Figure  42.  Convergence  of  the  LMS  algorithm  at  UAV  speed  of  80  m/s. 
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240  samples/s 
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Figure  43.  Mean  square  error  of  the  LMS  algorithm  at  different  UAV  speed. 
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H.  SUMMARY 

This  chapter  presented  the  simulation  model  and  the  results.  Simulation  results 
showing  the  effects  of  the  sensor  node  density,  position  errors  and  sensor  node  failures  on 
the  antenna  beam  generated  by  the  sensor  cluster  were  illustrated.  Results  showed  that 
the  average  sidelobe  levels  increased  while  the  maximum  average  power  gain  decreased 
when  the  density  (i.e.,  number  of  sensor  nodes  per  unit  area)  of  the  sensor  nodes  de¬ 
creased.  However,  the  shape  of  the  main  lobe  remained  relatively  unchanged  when  the 
density  changed.  Results  also  showed  that  the  presence  of  position  errors  and  sensor 
node  failures  reduced  the  maximum  average  power  gain  of  the  antenna  beam  and  in¬ 
creased  its  sidelobe  levels,  but  the  shape  of  the  main  lobe  remained  relatively  unchanged. 
These  simulation  results  were  found  to  be  in  good  agreement  with  the  theoretical  results 
for  random  arrays  [7,  8]. 

A  series  of  plots  showing  the  adaptive  beamforming  capability  of  the  sensor  clus¬ 
ter  were  also  illustrated.  It  was  shown  that  the  maximum  average  power  gain  and  the  in¬ 
terference  rejection  levels  decreased  when  the  position  errors  and  the  number  of  sensor 
node  failures  increased.  Additionally,  main  lobe  pointing  error  occurred  when  the  posi¬ 
tion  errors  increased  to  0.47,  and  the  number  of  sensor  node  failures  increased  to  twenty. 

Lastly,  the  significance  of  the  sampling  rate  on  the  convergence  of  the  LMS  algo¬ 
rithm  when  the  UAV  speed  increased  was  discussed. 
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V.  CONCLUSIONS 


A  wireless  senor  network  (WSN)  consists  of  a  large  number  of  small  sensor  nodes 
that  are  densely  deployed  over  an  area  of  interest  to  acquire  information  about  targets  of 
interest.  These  sensor  nodes  collaborate  among  themselves  to  form  an  ad-hoc  network 
and  disseminate  the  collected  target  information  to  an  UAV.  The  objective  is  to  increase 
the  data  rate  and  transmission  range  between  the  sensor  nodes  and  the  UAV.  To  meet 
this  objective,  a  distributed  beamforming  approach  was  proposed  whereby  the  sensor 
nodes  are  grouped  into  clusters  and  their  transmission  are  coordinated  in  order  to  form  a 
distributed  antenna  array  that  directs  a  beam  towards  the  UAV. 

A  three-tiered  hierarchical  clustering  sensor  network  architecture  was  considered 
in  this  thesis.  A  simulation  model  was  developed  and  implemented  in  MATLAB  pro¬ 
gramming  language  to  study  the  application  of  beamforming  using  sensor  clusters  to  es¬ 
tablish  a  communication  link  with  the  UAV.  Different  sensor  node  densities,  position 
errors  and  sensor  node  failures  were  simulated  to  investigate  their  effects  on  the  antenna 
beam  generated  by  the  sensor  cluster. 

A.  SIGNIFICANT  RESULTS 

Results  showed  that  the  average  sidelobe  levels  increased  while  the  maximum  av¬ 
erage  power  gain  decreased  when  the  density  (i.e.,  number  of  sensor  nodes  per  unit  area) 
of  the  sensor  nodes  decreased.  However,  the  shape  of  the  main  lobe  remained  relatively 
unchanged  when  the  density  changed.  Results  also  showed  that  the  presence  of  position 
errors  and  sensor  node  failures  reduced  the  maximum  average  power  gain  of  the  antenna 
beam  and  increased  its  sidelobe  levels,  but  the  shape  of  the  main  lobe  remained  relatively 
unchanged.  These  simulation  results  were  found  to  be  in  good  agreement  with  the  theo¬ 
retical  findings  for  random  arrays  [7,  8]. 

Since  the  shape  of  the  antenna  main  lobe  remains  relatively  unchanged  in  the 
presence  of  position  errors  and  sensor  node  failures  or  when  the  density  of  the  sensor 
nodes  changed,  it  can  be  steered  towards  the  UAV  in  an  adaptive  manner  as  the  UAV 
flew  over  the  sensor  field.  This  implied  that  the  energy  of  the  antenna  beam  can  be  di- 
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rected  towards  the  UAV,  resulting  in  higher  signal-to-noise  ratio,  and,  hence  higher  chan¬ 
nel  capacity.  Additionally,  the  maximum  average  power  gain  of  the  main  lobe  can  be  in¬ 
creased  by  increasing  the  density  of  the  sensor  cluster,  thereby  increasing  the  transmis¬ 
sion  range  between  the  sensor  clusters  and  the  UAV. 

B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

It  was  assumed  in  this  study  that  there  was  perfect  frequency,  phase,  data  syn¬ 
chronization  and  beamforming  coordination  among  the  sensor  nodes  within  the  cluster. 
This  may  be  difficult  to  achieve  in  practice.  The  effects  of  these  imperfections  may  be 
investigated  in  a  future  study. 

The  effects  of  mutual  coupling  on  the  antenna  beam  was  also  not  considered  in 
this  study.  In  reality,  some  sensor  nodes  may  be  close  to  each  other  and  the  effects  of 
mutual  coupling  on  the  antenna  beam  can  be  significant.  These  effects  on  the  quality  of 
the  antenna  beam  may  be  investigated. 

The  beamforming  process  was  implemented  on  the  sensor  clusters  in  this  study. 
The  same  idea  can  be  extended  to  the  UAV  so  that  the  performance  of  the  entire  MIMO 
system  can  be  evaluated. 


72 


APPENDIX 


MATLAB  PROGRAM  FLOWCHART 


This  Appendix  describes  the  flowchart  of  the  MATLAB  program  used  in  this  the¬ 
sis  work.  Figure  A.  1  shows  the  main  program  flowchart  and  the  sub-programs  are  de¬ 
scribed  in  the  following  discussions. 


Figure  A.  1 .  Main  program  flowchart 


Uniform  _1D. m 

This  program  simulates  the  normalized  power  gain  of  a  uniformly  excited  10x1  omni¬ 
directional  antenna  array  with  equal  element  to  element  spacing. 

Uniform  _2 DM  AIN  m 

This  program  simulates  the  normalized  power  gain  of  a  uniformly  excited  MxN  omni¬ 
directional  antenna  array  with  equal  element  to  element  spacing.  The  major  limitation  of 
this  program  is  that  the  M  must  be  equal  to  N. 

Random  _M A  IN.  m 

This  program  simulates  the  average  power  gain  of  a  random  array  with  MxN  elements 
subjected  to  random  positioning,  position  errors  and  sensor  node  failures.  The  major 
limitation  of  this  program  is  that  the  M  must  be  equal  to  N. 
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PG _plot_all.m 

This  program  produces  all  the  antenna  beam  pattern  plots  contained  in  this  thesis. 

Figures  A. 2  and  A. 3  show  the  flowchart  of  the  Uniform_2DMAIN.m  program  and 
the  and  the  Random  _M AIN. m  program,  respectively.  The  sub-programs  are  described  as 
followed. 

Un  iform_2DLOC.  m 

This  program  generates  the  x  and  y-coordinates  of  the  sensor  nodes  with  equal  element  to 
element  spacing. 

Uniform _2D  Wts.  m 

This  program  generates  the  magnitude  (all  equal)  and  angle  of  the  complex  weights. 
Uniform _2 DA  F_cal.  m 

This  program  computes  the  array  factor  of  the  sensor  cluster. 

PG _plot.m 

This  program  produces  the  antenna  beam  pattern  plot. 

Random  LOC. m 

This  program  generates  the  x  and  y-coordinates  of  the  sensor  nodes  with  random  posi¬ 
tions  within  the  specified  square  area  of  length  L. 

LMS.m 

This  program  computes  the  complex  weights  of  the  sensor  cluster  using  LMS  algorithm. 
Random  FA  IL.  m 

This  program  chooses  fail  number  of  elements  from  NxM  elements  and  set  the  magnitude 
of  their  complex  weights  to  zero. 

AF_cal.m 

This  program  computes  the  array  factor  of  the  sensor  cluster. 
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M  N.  wavelength  spacing,  theta 0.  phiO 


Figure  A. 2.  Program  flowchart  for  Uniform _2DM AIN. m 
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Na,  M,  N.  wavelength,  L,  error,  fail,  fs,  thetaO 


PG_plot.m 


Figure  A. 3.  Program  flowchart  for  Random _M AIN. m 
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PG _plot.m 

This  program  produces  the  antenna  beam  pattern  plot. 

Description  of  variables 

1 .  M=  number  of  elements  along  x-axis 

2.  N=  number  of  elements  along  y-axis 

3.  x  =  x-coordinates  of  the  sensor  node 

4.  xerr  =  x-coordinates  of  sensor  node  with  errors  (in  meters) 

5.  y  —  y-coordinates  of  the  sensor  node 

6.  y_err  =  y-coordinates  of  sensor  node  with  errors  (in  meters) 

7.  wavelength  =  design  wavelength  of  the  array  (in  meters) 

8.  L  =  length  of  square  plane  (in  meters) 

9.  error  =  position  errors  (in  meters) 

1 0.  fail  =  number  of  elements  failure 

1 1 .  spacing  =  element-element  spacing  along  x-axis  and  y-axis 

12.  thetaO  =  AOA  elevation  angle  (in  radians) 

13.  phiO  =  AOA  azimuth  angle  (in  radians) 

14.  Na  =  number  of  random  array  samples  to  be  generated 

15.  fs  =  sampling  rate 

1 6.  mag_Wts  =  magnitude  of  the  complex  weights 

1 7.  mag_Wts Jail  =  weights  with  magnitude  of  failed  elements  set  to  zero 

1 8.  ang_Wts  =  angle  of  the  complex  weights  (in  radians) 

1 9.  AFjnag  =  array  factor 

20.  avg  AF  mag  =  average  array  factor 
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